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1  Executive  Summary 

This  progress  report  summarizes  the  work  done  under  the  Grant  No.  N00014-89-K-2030 
entitled  "Experimental  and  Theoretical  Studies  of  High  Power  Plasma  Filled  Backward  Wave 
Oscillators"  during  the  period  August  1.  19S9  to  January  31,  1991. 

Relativistic  Backward  Wave  Oscillators  (BW'Os)  have  proven  to  be  relatively  efficient, 
high  power  microwave  sources  at  the  centimeter  and  millimeter  wavelength  range.  It  consists 
of  a  relativistic  electron  beam,  a  strong  guiding  magnetic  field  and  a  slow  wave  structure.  Our 
approach  for  increasing  the  power  generating  capabilities  of  relativistic  BWO’s  is  to  introduce 
plasma  into  the  device.  It  is  expected  that  in  this  way  very  large  electron  beam  currents 
(well  above  the  vacuum  space  charge  limit)  may  be  propagated,  resulting  in  extremely  high 
power  microwave  generation.  Also,  the  plasma  is  expected  to  reduce  the  ac  space  charge  in 
the  beam,  resulting  in  enhanced  axial  bunching  and  possibly  enhancement  in  efficiency. 

Previously,  it  was  found  experimentally  that  peak  efficiency  in  the  presence  of  the  plasma 
can  be  greatly  enhanced,  by  a  factor  of  8  (to  about  40%).  Also,  we  found  that  the  plasma  can 
be  externally  generated  and  then  injected  into  the  device.  The  experimental  configuration 
is  shown  in  Fig.  1.  For  a  detailed  drawing  of  the  slow  wave  structure  see  Fig.  2. 

The  purpose  of  the  research  covered  during  this  period  was  to  study  the  behavior  of 
a  plasma  loaded  BWO  when  driven  by  large  beam  current  approaching  the  vacuum  space 
charge  limiting  value  inside  the  device.  The  main  results  are: 

1.  The  plasma  indeed  allows  injection  of  large  current  into  the  BWO.  Currents  as  high 
as  75%  of  the  vacuum  space  charge  limit  were  successfully  propagated. 

2.  In  plasma  loaded  BWO,  in  contrast  to  the  vacuum  case,  the  enhanced  efficiency  can 
be  maintained  even  for  large  beam  currents  (approaching  the  vacuum  limit).  It  is 
anticipated  that  this  might  hold  even  well  beyond  the  vacuum  limit. 

3.  The  output  frequency  is  only  slightly  shifted  ( 1  %— o%)  bv  the  presence  of  the  plasma. 

4.  A  foilless  field  emission  diode  can  be  operate  successfully  in  a  background  of  low  density 
plasma. 

5.  An  important  parameter,  the  start  oscillation  current  was  measured  for  the  first  time 
in  vacuum  BWO  for  various  beam  voltages  and  geometries.  For  example,  at  600  kV 
the  start  current  is  under  140  A. 

6.  Initial  indications  are  that  the  plasma  increases  the  microwave  pulse  duration. 

7.  A  novel  field  emission  cathode  design  (new  geometries  and  materials)  allows  for  contin¬ 
uous  control  of  the  beam  current  in  the  range  of  50-6000  Amps.  Turn  on  fields  below 
•50  kY/cm  were  achieved. 

S.  Linear  theory  for  a  plasma  loaded  BWO  was  developed  (using  a  simplified  model). 

A  nonlinear  theory  for  a  finite  length  device  was  developed. 

We  achieved  considerable  progress  on  the  experimental  and  theoretical  aspects  of  microwave 
plasma  electronics,  as  summarized  below  and  reported  in  greater  detail  in  Section  2  of  this 
report. 


1:  Schematic  diagram  of  the  plasma  loaded  BVVO  experiment. 
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Figure  2:  The  rippled-wall  slow-wave  structure  used  in  our  i  asma  loaded  BVVO. 


Novel  Method  for  Determining  the  Electromagnetic  Dispersion  Relation  of  Pe¬ 
riodic  Slow  Wave  Structures  A  novel  method  for  calculating  the  dispersion  relation  of 
electromagnetic  modes  in  arbitrary  periodic  slow  wave  structures  was  found.  In  this  method 
it  is  sufficient  to  know  the  frequencies  corresponding  to  three  special  wave  number  values, 
with  other  points  calculated  using  an  approximate  analytical  expression.  This  technique  was 
successfully  applied  to  determine  the  dispersion  relation  of  T A/qi  mode  in  a  variety  of  slow 
wave  structures.  The  results  are  in  excellent  agreement  (within  0.15%)  with  other  numerical 
and  analytical  techniques.  (See  Section  2.1.) 

Studies  of  Novel  Field  Emission  Guns  for  BWOs  Experimental  studies  of  novel  non- 
planar,  magnetically  insulated,  foilless,  explosive  emission  guns  were  successfully  completed. 
New  geometries  and  materials  were  evaluated,  resulting  in  the  development  of  guns  which 
turn  on  at  very  low  average  electric  fields  (about  50  kV/cm).  These  geometries  inctuue 
annular,  single  cathode  design  and  annular,  cylindrically  symmetric  double  cathode  design. 
Additional  studies  investigated  the  effects  of  several  cathode  emitting  surfaces.  We  can  now 
vary  the  beam  current  in  the  range  of  50  A  to  6  kA  reproduciblv.  (See  Section  2.2.) 

The  Physics  of  Relativistic  Electron  Beams  Close  to  the  Vacuum  Space  Charge 
Limiting  Current  This  topic  is  important  in  the  understanding  of  intense  beam  interac¬ 
tion  in  BWOs.  It  was  found  that  the  potential  depression  due  to  dc  space  charge  drasti  ally 
change  the  behavior  of  the  space  charge  waves  supported  by  the  beams.  (These  waves  are 
responsible  for  the  interaction  in  a  BWO.)  For  example,  the  phase  velocity  of  the  slow  space 
charge  wave  vanishes  at  the  limiting  current.  It  was  also  found  that  the  beam  can  support 
four  space  charge  waves.  A  complete,  fully  relativistic  dispersion  relation  for  a  thin,  annular 
beam  in  an  infinitely  large  magnetic  field  was  analyzed  numerically.  (See  Section  2.3.) 

Plasma  Characterization  Measurements  were  made  of  plasma  density,  drift  velocity 
and  temperature  as  a  function  of  time,  space  and  magnetic  field.  Langmuir  orobes  and 
interferometric  techniques  were  used  to  characterize  two  kinds  of  plasma  gim«  -'ee  Section 
2.4.) 

Experimental  Studies  of  the  Starting  Current  in  a  Relativistic  Vacuum  BWO 

The  vacuum  BWO  was  studied  experimentally  over  wide  range  of  beam  currents,  from  the 
threshold  (starting  current)  to  75%  of  the  limiting  current  and  over  wide  range  of  voltages 
and  magnetic  fields.  Frequency,  power  and  pulse  shape  were  analyzed.  For  a  600  kV  beam, 
the  experimentally  measured  start  current  is  under  140  A.  (See  Section  2.5.) 

Frequency  Pulling  in  Relativistic  Vacuum  BWOs  The  frequency  of  a  vacuum  BWO 
can  be  tuned  bv  varying  the  beam  voltage.  Since  the  beam's  space  charge  waves  are  dras¬ 
tically  modified  by  space  charge  effects  close  to  the  limiting  current,  the  exact  relativistic 
dispersion  relation  was  used  to  calculate  the  new  interaction  frequency.  This  frequency  is 
defined  here  as  the  intersection  of  the  slow  space  charge  beam  wave  with  the  cold  structure 
dispersion  relation.  The  difference  between  the  two  is  called  "frequency  pulling".  It  was 
calculated  numerically  for  the  cases  of  interest.  (See  Section  2.6.) 


Experimental  Studies  of  Plasma  Loaded  BWOs  A  plasma  filled  relativistic  BWO  was 
driven  bv  an  intense  relativistic  electron  beam  with  current  approaching  the  vacuum  space 
charge  limiting  value.  The  beam  propagation  in  plasma  as  well  as  the  RVVO  characteristics 
"iich  as  power,  frequency,  and  pulse  duration  were  measured.  Initial  results  indicate  that 
the  presence  of  plasma  causes  small  frequency  shifts  (compared  with  the  vacuum  easel  in 
addition  to  efficiency  enhancement.  It  also  helps  to  extend  the  pulse  duration.  Above  a 
critical  plasma  density  the  interaction  is  quenched.  (See  Section  2.7.) 

Finite  Length  Effects  in  BWOs  When  finite  length  effects  (such  as  reflection  coefficients 
at  the  ends)  are  taken  into  accounts,  the  slow  wave  structure  behave  much  like  a  cavity  with 
discrete  axial  resonances  for  a  fixed  transverse  mode.  The  number  of  these  resonances  is 
equal  to  the  number  of  periods  in  the  slow  wave  structure  plus  one.  We  have  designed  and 
fabricated  a  cold  test  system  to  measure  these  resonances  as  well  as  the  reflection  coefficient 
and  the  dispersion  relation  of  the  TM0i  mode  in  the  actual  structure.  These  resonances  are 
due  to  various  axial  modes.  (See  Section  2.8.) 

Linear  Theory  of  Plasma  Loaded  BWO  Linear  theory  of  plasma  loaded  BWO  was 
developed  (for  infinitely  strong  guiding  magnetic  field,  sofid  beam  and  neglecting  electron 
plasma  waves).  It  shows  that  below  some  critical  background  plasma  density  the  instability 
is  absolute  (BWO),  while  above  it  is  convective  (TWT)  in  nature.  Growth  rates  were  also 
calculated.  The  linear  theory  will  have  to  be  expanded  to  better  explain  the  experimental 
results.  (See  Section  2.9.) 

Theory  of  Relativistic  BWOs  As  a  first  step  in  the  development  of  a  comprehensive 
theory  of  plasma  loaded  relativistic  BWO,  we  have  developed  a  linear  and  nonlinear  theory 
for  a  vacuum  relativistic  BWO  of  a  finite  length.  The  main  result  of  our  investigations  is 
the  discovery  of  regions  in  parameter  space  of  single  and  multimode  BWO  operation.  We 
found  that  the  start  oscillation  current  and  the  efficiency  of  a  single  frequency  operation  are 
sensitive  to  the  choice  of  a  operation  point  of  a  BWO.  The  operation  point  in  stationary 
regime  is  modified  by  the  presence  of  a  low  density  plasma.  However,  preliminary  estimates 
indicate  that  Maryland  BWO  is  operating  in  overbunch  unstable  region  (multimode  regime). 

I  lie  eject  of  plasma  in  this  regime  will  be  a  subject  of  the  future  investigations.  (See  Section 

2.10.) 

Open  Research  Topics  We  achieved  considerable  progress  on  both  the  experimental  and 
theoretical  aspects  of  plasma  loaded  BWO.  Additional  research  is  needed,  however  in  order 
to  take  full  benefit  of  the  innovative  approach  of  plasma  microwave  electronics.  In  particular 
we  need  to  address  such  questions  as  operation  well  and  above  the  vacuum  space  charge 
limiting  current,  over  moded  systems  and  understanding  the  operating  regimes  and  pulse 
duration.  iSee  Section  2.11.) 

Our  results  were  extensively  documented,  published  in  a  timely  manner,  and  caused 
<  onsiderabie  scientific  interest.  Fifteen  papers  of  various  kinds  were  published,  including  two 
in  Physical  Review  Letters  and  three  invited  talks.  (See  Section  2.18.) 


2  Research  Progress 

2.1  Novel  Method  for  Determining  the  Electromagnetic  Dispersion  Relation  of 
Periodic  Slow  Wave  Structures 

It  is  important  to  accurately  know  the  dispersion  relation  of  the  slow  wave  structure  in  order 
to  synchronize  the  phase  velocity  of  the  electromagnetic  wave  with  the  electron  beam  to 
produce  efficient  interaction.  There  are  few  techniques  to  calculate  the  dispersion  relation  of 
a  slow  wave  structure.  One  technique  relies  on  an  equivalent  lumped-element  circuit  model. 
Others  are  analytical  and  are  suitable  for  sinusoidally  or  square  wall  small  perturbations. 
Numerical  techniques  (like  the  code  SUPERFISH)  require  considerable  effort  and  time  in 
order  to  accurately  calculate  the  dispersion  relation  everywhere  in  the  u-k  plane.  Here,  we 
developed  a  novel,  accurate  and  simple  method  for  calculating  the  dispersion  relation  of  a 
slow  wave  structure  of  arbitrary  geometry. 

•  Basic  Principle.  For  any  periodic  slow  wave  structure  with  n  periods,  when  shorted 
at  both  ends,  there  are  (n  4-  1)  discrete  resonant  frequencies  with  a  phase  shift  per 
period  equally  spaced  between  0  and  jr.  Instead  of  using  a  large  number  of  points  to 
accurately  calculate  the  periodic  dispersion  relation  in  the  range  0  to  we  calculate 
the  resonant  frequencies  corresponding  to  three  special  normalized  wavenumber  values 
and  then  use  an  analytically  expression  to  produce  the  complete  dispersion  curve. 

•  Results.  For  our  sinusoidally  corrugated  slow  wave  structure  (see  Fig.  2);  the  radius  is 
described  by  R  =  Ro[l  +  h  cos(2 icz/L)\  where  Rq  is  the  average  radius,  h  is  the  normal¬ 
ized  corrugation  amplitude,  and  L  is  the  structure  period.  Superfish  modeling  yielded 
(with  proper  boundary  conditions)  the  three  special  resonance  frequencies  /0,  /*•/ 2,  /* 
in  the  lower  passband.  The  results  are  tabulate  in  Table  1.  The  corresponding  mode 
patterns  (electric  field  lines)  are  shown  in  Fig.  3. 

The  approximate  form  of  the  dispersion  relation  is  given  by  /  =  A  —  B  cos  3L  —  C  cos*  3L. 
where  A  =  B  =  (/,  -  f0)/2:C  =  fx/2  -  (/*  +  fo)/~-  The  calculated  dispersion  relation 
is  plotted  in  Fig.  4. 

In  order  to  estimate  the  accuracy  of  this  technique,  we  first  calculate  additional  points 
on  the  dispersion  curve  using  SUPERFISH.  As  an  independent  check,  the  dispersion  relation 
derived  here  was  compared  to  that  derived  by  other  techniques  for  the  same  structure.  Both 
are  in  excellent  agreement  (within  0.15%)  over  the  entire  frequency  band. 


FREQUENCY  (GHz) 


able  1:  The  resonance  frequencies  (GHz)  of  a  sinusoidally  corrugated  waveguide  cavitv 

(Ro  =  1.5  cm,  L  =  1.67  cm.  h  =  0.273). 
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Comparison  of  the  TM0i  dispersion  curve  of  other  technique  (solid  line)  and  this 
woik  (dotted  line)  for  a  sinusoidally  corrugated  waveguide. 


2.2  Studies  of  Novel  Field  Emission  Guns  for  BWOs 

Implosive  emission  diodes  are  often  used  to  generate  the  intense  relativistic  electron  beams 
used  in  high  power  microwave  devices.  In  contrast  to  typical  explosive  emission  nuns  which 
produce ■  beams  with  current  densities  on  the  order  10  k A/cnr.  we  have  studied  cathode 
designs  m  which  the  beam  current  density  is  tunable  over  a  wide  range  ('200  A/cnr  -12.000 
A/cnr  i.  This  is  necessary  in  order  to  access  different  operating  regime  of  BWO.  from  the 
threshold  (starting  current)  and  up. 


Purpose  of  This  Work 

•  1’se  magnetically  insulated,  foilless  emission  explosive  diodes  to  generate  a  hollow 
relativistic  electron  beam  in  a  BWO  slow  wave  structure.  The  beam  radii  are  7.75 
mm  and  2.6  mm  with  thicknesses  0.1  to  0.2  cm. 

•  Carefully  analyze  the  operation  of  explosive  emission  diodes  to  determine  how  the 
microwave  output  power  and  resonant  frequency  depend  on  the  parameters  of  the 
beam. 

•  Study  the  characteristics  of  the  gun  injection  and  beam  propagation  in  magnetized 
BWO  structures. 

•  Scan  a  wide  range  of  currents  from  the  BWO  start  oscillation  current  to  the  space 
charge  limiting  current.  Both  single  and  double  cathode  geometries  are  employed  (see 
Fig.  5). 

•  Analyze  the  beam  uniformity. 

•  Study  cathode  coating  techniques  to  reduce  the  average  electric  field  required  for  emis¬ 
sion. 


The  Relativistic  Diode  Analysis  Suppose  that  the  anode  potential  is  zero,  while  t lie 
•  uihode  is  at  potential  of  On.  Conservation  of  energy  implies  that 


=  1  -t- 


eo 

mnc: 


By  solving  Poisson  s  equation  with  the  boundary  condition: 


,1-0  .! 

fol,o 

We  get  the  following  relationship  between  current  density  and  the  applied  voltage  for  the 
relativistic  planar  diode: 


rnnc ft, 
'led1 
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■  iMw  /■  anti  E  are  elliptic  integrals  of  the  first  and  second  kind.  In  the  limit  of  nonreiativistic 
■earn.  io/m0c:  <<  1.  Eq.|3)  reduces  to  the  familiar  ( 'hild-Langmuir  expression: 

t  '  e.  \  1/2  o3/2 

=  ij  t  mo  J  qr 

In  the  ultra-relativistic  limit.  *  2>  1.  Eq.(3)  reduces  to 

rr2  =  i/2  'hi 

where  /'  =  fo/mnc2  and  C*  —  <  d:J/m0c^Co. 

A  comparison  of  the  normalized  diode  impedance  for  the  relativistic  and  non-relativistic 
-elution  is  presented  in  Fig.  6.  The  diode  impedance  decreases  as  the  applied  potential  is 
increased.  And  also,  the  calculated  current  as  a  function  of  diode  voltage  is  shown  in  Fig.  '> 
'or  a  certain  AK  gap  1 12.5  mm). 

Conclusions 

1.  The  gun  followed  a  Child-Langmuir-like  law  (I  ~  V'or).  The  experimentally  measured 
power  coefficient  varies  from  1  to  1.5  for  different  gun  geometries,  compared  with  1.5 
for  the  non-relativistic  Child  Langmuir  law.  Figures  7  and  8  show  the  experimental 
V-I  characteristic  of  different  gun  geometries. 

2.  t’sing  both  single  and  double  cathode  configurations  (see  Fig.  5).  the  beam  current 
ranges  from  100  A  to  0  k.A  (current  density  from  200  A/cm2  to  12  kA/cnr  ).  fleam 
•  ■nergies  from  200  keV  to  800  keV  were  generated. 

3.  The  effective  emission  area  is  larger  than  the  physical  cathode  area.  For  example,  fur 
AK  Cap=12.5  mm.  .  l„(r  =  10.8.-lphvs  i  see  Fig.  0). 

1.  ( Iraphite/epoxy  cathode  coating  reduced  the  average  electric  field  required  for  uniform 
'•mission  to  about  50  kV/crn.  With  the  coating,  we  obtain  a  uniform  beam  at  150  A. 
while  without  the  coating,  this  value  was  350  A  (400  kV)  and  the  beam  was  filamented. 


Figure  5:  Schematic  diagram  of  the  gun  geometries  as  well  as  representative  damage 

patterns 
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Norllzftd  Impedance  (Zn=Uf  ) 


Planar  Diode  Gun  Calculation 


Norlized  Voltage  (U=eV/(mc') 


F igure  6:  A  comparison  of  the  relativistic  and  nonreiativistic  solution  for  the  impedance  of 

a  one-dimensional  planar  diode. 
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2.3  The  Physics  of  Relativistic  Electron  Beams  Close  to  the  Vacuum  Space 
Charge  Limiting  Current 


I  his  subject  is  important  for  the  understanding  os  intense  beams  interaction  in  BWOs. 
especially  when  the  beam  current  approach  the  limiting  value.  This  is  clearly  the  case  here, 
since  under  these  condition  we  might  expect  to  maximize  the  microwave  output  power  (the 
plasma,  of  course,  is  expected  to  drastically  change  the  beam  behavior). 

The  electron  beam  drift  energy  is  influenced  by  the  diode  potential  (7 c),  beam  current 
density  ( L,)  and  the  beam  geometry.  By  solving  Poisson's  equation  with  the  appropriate 
boundary  condition,  we  get  the  well  known  dependence  of  /beam  on  the  effective  beam  energy 
•:b  is  as  follows 


h 


U 

2ln(r0/rb) 


(7i-l)l/2(7c/7b-U 


(') 


where  r0  is  the  drift  tube  radius  and  77  is  the  beam  radius.  Maximizing  this  expression 
with  respect  to  the  injection  energy  yield  the  absolute  maximum  current  which  can  stably 
propagate  under  the  ideal  condition  of  infinitely  large  guiding  magnetic  field.  This  maximum 
value  is  plotted  in  Fig.  11,  while  the  beam  geometry  is  defined  in  Fig.  10. 

\\  hen  an  electron  beam  is  injected  into  the  slow  wave  structure  (instead  of  a  smooth  wall 
waveguide/drift  tube)  it  may  excite  an  unstable  interaction.  This  instability  occurs  when 
a  normal  mode  of  structure  having  positive  energy,  interacts  with  the  negative  energy,  slow 
space-charge  wave  of  electron  beam.  By  varying  the  slope  of  the  slow  space-charge  wave  of 
the  electron  beam,  the  interaction  point  is  shifted  along  the  dispersion  curve  of  the  structure. 
To  accurately  predict  the  operating  frequency  of  the  BWO,  an  understanding  of  dispersive 
characteristics  of  the  slow  space  wave  of  the  beam  close  to  the  vacuum  space  charge  limiting 
current  is  necessary. 

The  dispersion  relation  for  the  space  charge  waves  derived  by  the  linear  wave  analysis  for 
an  annular,  infinitely  thin  relativistic  electron  beam  could  be  expressed  as  follows: 


■/o(rr0)  4-  a 


Tc 

ui  —  kvb 


(rrb){J0(rrb)N0{Tr0)  -  J0[Tr0)N0(Trb)]  =  0 


(8) 


where  a  =  r /#,/ ^7^ / 4 ,  T2  =  ■S' / c2  —  k2  and  -V0  are  the  zeroth  order  Bessel  functions  of 
the  first  and  second  kind,  respectively,  a  represents  the  space  charge  parameter. 

Solving  this  equation  numerically,  we  could  get  the  u  vs  k  curves  for  the  two  space  charge 
wave,  the  fast  (positive  energy)  and  the  slow  (negative  energy).  In  Fig.  12.  we  show  only  the 
-low  space  charge  wave  (labeled  SSCW)  for  various  beam  currents.  We  clearly  see  from  this 
figure  that  as  the  current  is  increased  (a  is  larger),  the  phase  velocity  u ;/k  decreases  for  the 
-tune  injection  energy.  This  phase  velocity  (w’/T'U— 0)  is  plotted  versus  the  beam  current  and 
it  vanishes  when  the  beam  current  exactly  equals  the  space  charge  limiting  value  in  vacuum 
1  see  Fig.  13).  Furthermore,  this  phase  velocity  can  also  be  negative.  To  understand  this 
-eemingly  strange  phenomena  let's  look  at  Fig.  14a,  which  is  a  plot  of  />,  vs  ~b-  For  a  given 
cathode  potential  and  a  certain  current,  7^  has  two  equilibrium  solutions.  One  corresponds 
to  a  higher  density  and  a  lower  velocity.  Another  corresponds  to  a  lower  density  and  high 
velocity.  These  two  equilibrium  solutions  merges  when  the  beam  current  exactly  equals  the 
vacuum  space  charge  limiting  value  (i.e..  the  peak  in  lb).  At  this  point,  the  beam  drift  energy 


I'b)  equals  to 


lb  =  7  c1/3  (9) 

We  see  that  there  are  two  equilibrium  solution  for  a  beam  in  a  smooth  drift  section  (i.e. 
waveguide),  and  therefore  the  beam  can  support  two  pairs  of  space  charge  waves.  As  an 
example,  lets  look  at  the  4  waves  supported  by  a  2  kA  beam  with  an  injection  energy  of 
600  kV.  The  conventional  fast  and  slow  wave  are  recovered  for  ~/f)  >  7^3  (Fig.  15)  and  a 
second  pair  exists  for  7 4  <  7C1/3  (Fig.  16).  The  solution  which  correspond  to  higher  beam 
density  and  lower  velocity  is  usually  labeled  as  unstable.  In  this  region,  the  electrostatic 
potential  energy  of  each  electron  is  larger  than  its  kinetic  energy.  As  indicated  earlier,  when 
the  beam  current  exactly  equal  to  the  space  charge  limiting  current,  the  two  solutions  merge 
(see  Fig.  17).  The  consequences  of  this  analysis  will  be  further  discussed  in  Section  2.6. 


I 

I 
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Figure  10:  The  configuration  of  the  beam 


geometry  under  consideration. 


gure  11:  Space  charge  limiting  current  vs  diode  voltage  with  drift  tube  radius  of  1.5  cm. 
beam  radius  of  0.775  cm  and  0.26  cm  (annular  beam). 
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[  inure  i 4:  The  reiation  between  the  beam  current,  the  space  charge  parameter  (a),  and  -\j. 
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Figure  15:  The  dispersion  relation  of  the  two  space  charge  waves  supported  by  a  beam  for 
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Figure  16:  The  dispersion  relation  of  the  two  space  charge  waves  supported  by  a  beam  tor 
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2.4  Plasma  Characterization 


In  order  to  understand  the  physical  mechanism  of  efficiency  enhancement  in  a  plasma  BWO. 
i .  is  important  to  accurately  measure  the  plasma  characteristics.  The  gun  is  located  in  a 
field  free  region  and  the  plasma  is  injected  along  the  axis  of  a  solenoidal  magnetic  field 
as  shown  in  Fig.  18a.  Measurements  were  made  of  the  plasma  density,  drift  velocity  and 
temperature  as  a  function  of  time,  space  and  magnetic  field.  We  have  examined  three 
types  of  plasma  guns  (Marshall.  Coaxial,  and  Hydrogen  Flashover  gun).  A  sketch  of  the 
geometry  for  two  of  these  two  guns  is  given  in  Fig.  18b, c.  We  used  both  Langmuir  probes 
and  microwave  interferometers  to  measure  the  plasma  density  with  and  without  magnetic 
field.  The  sketch  of  the  interferometer  block  diagram  and  Langmuir  probe  measurement 
circuit  diagram  was  shown  in  Fig.  19.  The  plasma  density  vs  distance  is  shown  in  Fig.  20. 
The  plasma  temperature,  plasma  drift  velocity  and  plasma  temperature  for  both  guns  are 
tabulated  in  Table  2  for  the  case  of  zero  magnetic  field.  With  magnetic  field,  the  Langmuir 
probe  was  located  about  1  m  away  from  the  gun  which  is  about  the  same  position  of  the 
interaction  region  of  the  BWO  experiment.  Figure  21  is  the  typical  result  of  plasma  density 
vs  time  inside  the  magnetic  field  (B  =  12  kG)  for  Hydrogen  Flashover  gun.  The  hydrogen 
:!ashover  gun  generated  a  plasma  cloud  moving  at  an  average  velocity  of  about  5-10  cm/ nsec 
and  with  a  temperature  of  5-15  eV.  The  detailed  results  are  tabulated  in  Table  3. 


SOLENOID  LANGMUIR  PROBE 


1  inure  18a:  A  schematic  diagram  of  the  plasma  gun  and  Langmuir  probe  position  relative 
to  the  solenoidal  magnetic  field  (for  plasma  characterization). 
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(b)  Marshall  Plasma  Gun 


(c)  Hydrogen  Plasma  Gun 


TiH2  in  Sodium  Silicate 
Solution 


Unit  in  inch 


^ure  18:  Sketch  of  geometry  for  (b)  Marshall  gun,  (c)  Hydrogen  Flashover  gun. 
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Figure  19:  (a)  Block  diagram  of  Interferometric  measurement,  (b)  Circuit  sketch  of 

Langmuir  Probe  measurement. 
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Figure  20:  Plasma  density  vs  distance  without  magnetic  field. 


Temperature  T  &  few  eV 

Plasma  drift  velocity  5  —10  cm/fisec 

1:  Tvpical  result  of  plasma  density  vs  time  inside  the  magnetic  Held  ( B  =  12  kG) 

for  Hydrogen  Flashover  Gun. 


Table  2:  Plasma  Characteristics  Measurement  Using 
Langmuir  Probe  Technique 


Hydrogen  Flashover  Plasma  Gun,  B  =  0 


Langmuir  Probe  located  41  cm  from  the  Gun 


Gun  Charge  Drift  Temperature  Plasma 

Voltage  Velocity  density 

(kV)  (cm//is)  (eV)  (cm-3) 


3.30  x  1012 
3.61  x  1012 
3.46  x  1012 


Langmuir  Probe  located  119  cm  from  the  Gun 


Gun  Charge 
Voltage 
(kV) 

Drift 
Velocity 
(cm/  fxs) 

Temperature 

(eV) 

Plasma 

Density 

(cm-3) 

10 

7.69 

3.33 

2.34  x  1010 

12 

8.33 

3.50 

3.74  x  10l° 

15 

8.50 

3.70 

5.40  x  10l° 

Marshall  Type  Plasma  Gun,  B  =  0 


Langmuir  Probe  located  119  cm  from  the  Gun 


Gas 

Type 

Gun  Charge 
Voltage 
(kV) 

Gas 

Pressure 

(psi) 

Temperature 

(eV) 

Plasma 

Density 

(cm-3) 

Argon 

5 

8 

3.85 

3.44  xlO 

Argon 

5 

10 

3.40  x  10 

Argon 

5 

12 

2.31  xlO 

Argon 

6 

8 

1.15  xlO 

Argon 

6 

10 

1.03  xlO 

Argon 

6 

12 

0.86  xlO 

Helium 

5 

8 

3.57 

1.43  xlO 

Helium 

5 

10 

0.95  xlO 

Helium 

5 

12 

1.07  xlO 

Nitrogen 

5 

8 

2.0 

2.55  xlO 

Nitrogen 

5 

10 

2.79  xlO 

Nitrogen 

5 

12 

2.71  xlO 

Table  3:  Plasma  Characteristics  Measurement  Using 
Langmuir  Probe  Technique 

Hydrogen  Flashover  Plasma  Gun 


Langmuir  Probe  located  119  cm  from  the  Gun 


Magnetic 

Field 

(kG) 

Gun  Charge 
Voltage 
(kV) 

Drift 
Velocity 
(cm/  ns) 

Temperature 

(eV) 

Plasma 

density 

(cm-3) 

1 

12.0 

10.6 

9.0 

2.70x  10u 

12.0 

11.2 

9.3 

2.40x10" 

12.0 

12.3 

10.2 

2.75x10" 

WEEm 

12.0 

12.5 

10.5 

2.69x10" 

Marshall  Type  Plasma  Gun 


Langmuir  Probe  located  119  cm  from  the  Gun 
Gun  Charge  Voltage  5  kV 


Gas  Type 

Magnetic 

Gas 

Temperature 

Plasma 

Field 

Pressure 

Density 

(kG) 

(psi) 

(eV) 

(cm-3) 

Argon 

9.0 

8 

5.0 

4.0x10" 

Argon 

10.4 

8 

4.0x10" 

Argon 

11.3 

8 

2.2x10" 

Argon 

12.25 

8 

5.4x10" 

Argon 

13.5 

8 

4.0x10" 

Argon 

15.5 

8 

4.0x10" 

Helium 

9.0 

8 

7.0 

2.3x10" 

Helium 

10.4 

8 

3.4x10" 

Helium 

11.3 

8 

3.5x10" 

Helium 

12.25 

8 

2.5x10" 

Nitrogen 

9.0 

8 

8.33 

4.3x10" 

Nitrogen 

10.4 

8 

3.4x10" 

Nitrogen 

11.3 

8 

4.3x10" 

Nitrogen 

12.25 

8 

4.3x10" 
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2.5  Experimental  Studies  of  the  Starting  Current  in  a  Relativistic  Vacuum 
BWO 

The  vacuum  BWO  was  studied  experimentally  over  wide  range  of  beam  currents,  from  the 
threshold  (starting  current)  to  75%  of  the  limiting  current  and  over  wide  range  of  voltages  and 
various  magnetic  fields.  Frequency,  power  and  pulse  shape  were  analyzed.  Most  important, 
however,  axe  the  measurements  of  the  BWO  starting  current.  Usually,  it  is  very  difficult 
to  achieve  low  beam  currents  with  field  emission  guns  which  favors  high  current  density 
(~10kA/cm2).  Using  the  novel  techniques  described  in  section  2.2,  we  were  able  to  reduce 
the  total  beam  current  to  a  value  below  100A  while  still  maintaining  a  good  beam  uniformity. 
This  enable  us  for  the  first  time  to  make  four  independent  measurements  of  the  starting 
current  (two  beam  diameter,  two  beam  energies).  The  experimental  results  are: 


1.  rb  =  0.775  cm 
I  start  <  114  A 

htart  <  HO  A 


at  430  kV  double  cathode  gun  AKI/O  =  53/50  mm 
at  600  kV  double  cathode  gun  AKI/O  =  58/55  mm 


2.  rb  =  0.26  cm 

114  A  <  I 3tart  <  335  A  at  430  kV 
132  A  <  I, tan  <  432  A  at  600  kV 


double  cathode  gun  AKI/O  =  55/55  mm 
double  cathode  gun  AKI/O  =  35/55  mm 


From  the  experiment,  it  was  found  that  as  the  beam  current  increases  above  /J(ort,  the 
envelope  of  the  microwave  pulse  is  effected  and  the  envelop  is  characterized  by  more  struc¬ 
tured.  At  /  >  I start  region,  we  found  even  more  structure  and  the  microwave  output  pulse 
is  narrower.  This  may  be  due  to  the  self-modulation  (see  section  2.10). 


28 


2.8  Frequency  Pulling  in  Relativistic  Vacuum  BWOs 

The  frequency  of  interaction  in  a  vacuum  BWO  is  expected  to  depend  on  the  beam  energy. 
Since  the  beam’s  slow  space  charge  wave  is  drastically  modified  by  space  charge  close  to 
the  limiting  current  (see  section  2.3)  the  exact  relativistic  dispersion  relation  was  used  to 
calculate  the  new  interaction  frequency.  This  frequency  is  defined  here  as  the  intersection 
of  the  SSCW  with  the  cold  TMoi  dispersion  relation  in  the  corrugated  waveguide.  The 
difference  between  the  two  is  called  “  frequency  pulling”.  Figure  22  shows  the  cold  slow 
wave  structure  dispersion  relation  superimposed  on  the  beam’s  slow  space  charge  wave  as 
calculated  for  various  beam  currents.  From  this  plot,  we  see  that  not  only  the  slope  of  the 
slow  space  charge  wave  is  changed  by  varying  the  beam  current  but  the  line  is  also  displaced. 
The  interaction  frequency  could  be  determined  by  the  intersection  point  of  the  dispersion 
curve  of  the  slow  space-charge  wave  and  the  TMot  structure  wave.  Figure  23  shows  the 
interaction  frequency  of  TMoi  mode  vs  the  diode  voltage  for  various  beam  currents.  For  the 
calculations  presented  in  this  figure  the  beam  radius  is  0.75  cm.  Figure  24  shows  similar 
calculation  for  a  beam  with  a  smaller  radius  (0.26  cm).  Figure  23  shows  that  the  interaction 
frequency  is  a  function  of  the  beam  current  even  for  a  fixed  gun  voltage,  i.e.  it  is  a  function 
of  the  gun  geometry  (like  anode-cathode  gap).  By  using  the  measured  relation  between  the 
gun  voltage  and  current  (see  section  2.2),  we  can  now  transform  Fig.  23  to  Fig.  25  which 
show  the  calculated  interaction  frequency  for  seven  different  experimental  gun  geometries 
which  were  studied  separately  (see  section  2.2).  The  horizontal  dashed  lines  correspond  to 
the  upper  and  lower  cut-off  frequencies  (band  pass)  of  the  TM0i  mode  in  the  slow  wave 
structure.  Figure  26  shows  the  experimentally  measured  frequencies  for  two  of  the  seven  gun 
geometries. 
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Figure  22:  Cold  slow  wave  structure  dispersion  relation  superimposed  on  the  beam’s  slow 
space  charge  wave  as  calculated  for  various  beam  currents. 
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igure  23:  The  interaction  frequency  of  TM0i  mode  vs  the  diode  voltage  for  various  beam 

currents,  n,  =  0.775  cm. 
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Figure  24:  The  interaction  frequency  of  TMoi  mode  vs  the  diode  voltage  for  various  beam 

currents,  rj  =  0.26  cm. 


Figure  25:  Calculated  interaction  frequencies  for  seven  different  experimental  gun 

geometries. 
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Figure  26:  The  experimentally  measured  frequencies  for  two  of  the  seven  gun  geometries. 


2.7  Experimental  Studies  of  Plasma  Loaded  BWOs 

A  plasma  filled  relativistic  BWO  was  driven  by  an  intense  relativistic  electron  beam  with 
current  approaching  the  vacuum  space  charge  limiting  value.  The  beam  propagation  in 
plasma  as  well  as  the  BWO  characteristics  like  power,  frequency  and  pulse  duration  were 
measured.  Initial  results  indicate  that  the  presence  of  plasma  causes  only  small  frequency 
shifts  (compared  with  the  vacuum  case)  in  addition  to  efficiency  enhancement.  It  also  helps 
to  extend  the  pulse  duration.  Above  a  critical  plasma  density  the  interaction  is  quenched. 

First  we  look  at  the  effect  of  plasma  on  beam  propagation.  For  anode-cathode  gap  of 
5.4  mm.  the  measured  gun  current  at  600  kV  is  about  4.5  kA  (beam  radius  is  0.75  cm.  see 
Fig.  5).  Under  the  same  conditions  the  calculated  space  charge  limiting  current  in  the  slow 
wave  structure  is  about  7.5  kA,  and  only  about  60%  of  the  injected  current  propagated.  As 
the  background  plasma  density  is  raised,  higher  and  higher  percentage  of  the  injected  current 
propagates.  At  a  background  plasma  density  of  about  2  x  1011  cm'3,  90%  of  the  beam  is 
transmitted  (see  Fig.  27).  An  enhancement  of  the  microwave  power  generation  efficiency 
was  observed  over  a  wide  range  of  injected  plasma  densities.  Variation  of  the  firing  delay 
between  the  plasma  gun  and  the  electron  beam  allowed  a  measure  of  control  over  the  plasma 
density  inside  the  slow  wave  structure.  The  microwave  peak  power  output  as  a  function  of 
this  firing  delay  is  shown  in  Fig.  28.  (These  results  were  achieved  with  the  Marshall  plasma 
gun.)  Microwave  output  rises  to  a  peak  of  almost  600  MW,  corresponding  to  an  electron 
efficiency  of  about  of  40%,  compared  with  about  5%  for  the  vacuum  BWO.  Under  similar 
conditions  this  enhanced  efficiency  is  maintained  in  a  plasma  loaded  BWO  even  for  large 
beam  currents  approaching  the  vacuum  space  charge  limit  (a  Hydrogen  Flashover  Gun  was 
used.) 

The  presence  of  the  plasma  caused  only  a  small  frequency  shift  compared  to  the  vacuum 
case.  This  frequency  shift  was  measured  by  heterodyning  the  BWOs  frequency  against  a 
preset  local  oscillator.  The  beat  frequency  between  the  two  signals  was  measured  on  a  high 
speed  recorder  (see  Fig.  29). 

Based  on  our  results,  we  can  draw  the  following  conclusions: 

1.  A  plasma  can  be  generated  externally  and  injected  into  the  BWO  interaction  area  in 
a  reproducible  and  controlled  manner. 

2.  A  plasma  filled  BWO  is  characterized  by  high  peak  efficiency  (about  40%)  and  small 
( 1  to  5%)  frequency  upshifts. 

1.  The  plasma  allows  for  the  injection  of  large  beam  currents.  In  our  experiments,  beam 
currents  of  up  to  75%  of  the  vacuum  limit  were  propagated.  Further  increases  well 
above  the  vacuum  limit  seem  feasible. 

1.  The  presence  of  the  plasma  increases  the  microwave  pulse  duration. 

It  appears  that  this  enhanced  efficiency  can  be  maintained  for  beam  currents  approaching 
and  exceeding  the  vacuum  limit.  A  slightly  overmoded  (diameter/wavelength  =  2  -  3)  will  be 
needed  for  peak  power  handling  capabilities  of  5  -  10  GW.  Further  research  may  concentrate 
on  investigating  such  systems. 

It  is  anticipated  that  plasma  injection  may  also  prove  beneficial  to  i  variety  of  high  power 
microwave  devices  such  as  gvrotrons  and  free  electron  lasers. 


Ib/Id  Beam  Current  /  Diode  Current 


Time 


CQ 

(U 


PQ 


Figure  29:  Beat  frequency  signal  vs 
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2.8  Finite  Length  Effects  in  Slow  Wave  Structures 

When  finite  length  effects  are  taken  into  account,  the  slow  wave  structure  behave  like  a  cavity 
with  discrete  resonances.  This  is  especially  pronounced  when  the  reflection  coefficients  at 
both  ends  are  large.  These  resonances  are  due  to  various  axi?l  modes.  The  number  of 
these  resonances  is  equal  to  the  number  of  periods  in  the  slow  wave  structure  plus  one. 
We  have  designed  and  fabricated  a  cold  test  system  to  measure  these  resonances  as  well 
as  the  reflection  coefficient  and  the  dispersion  relation  of  the  TM0i  mode  in  the  slow  wave 
structure.  The  cold  test  system  is  shown  in  Fig.  30,  and  will  be  excited  by  a  probe  (labeled 
6  in  the  figure).  The  structure  is  modular  and  can  vary  in  length  from  four  to  six  periods. 
The  expected  resonances  for  a  four  period  and  six  period  system  are  given  in  Table  4.  The 
results  will  be  compared  with  the  calculation  (see  Section  2.1). 
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Figure  30:  The  schematic  drawing  of  the  cold  test  system. 


Table  4:  Frequencies  Expected  from  Cold  Test 
of  the  Slow  Wave  Structure  (GHz) 


Sinusoidally  Corrugated  Waveguide 
( Ro  —  1.5  cm.  L  =  1.67  cm.  h  =  0.273) 


Resonance  Frequency 

4  Period  Cavity 

6  Period  Cavity  | 

fo  =  f(3L  =  0) 

7.404 

7.404 

f*,«  =  fiiSL  =  rr/6) 

N/A 

7.448 

!  fr/ 4  =  f(31‘  =  Z  /  A) 

7.558 

N/A 

f*l 3  =  f{0L  =  x/3) 

N/A 

7.702  : 

f.n  =  f(3L  =  z/  2) 

8.046 

8.046 

fij/z  =  /iTL  =  2-/3) 

N/A 

8.386 

I  h t/4  =  f(JL  =  371-/4) 

S.526 

N/A 

,  hn,e  =  f(3L  =  5?r/6) 

N/A 

8.632  , 

/*  =  flJL  =  T) 

8.773 

8.773 
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2,9  Linear  Theory  of  Plasma  Loaded  BWO 

A  linear  theory  of  plasma  loaded  BWO  was  developed  (for  infinitely  strong  guiding  magnetic 
field,  solid  beam  ).  It  shows  that  below  some  critical  background  plasma  density  the  insta¬ 
bility  is  absolute  (BWO),  while  above  it,  it  is  convective  (TWT)  in  nature.  Growth  rates 
were  also  calculated. 

The  linear  model  was  derived  for  the  excitation  of  electromagnetic  waves  in  a  plasma- 
filled  corrugated- wall  waveguide  with  an  arbitrarily  large  sinusoidal  corrugation.  The  model 
neglects  the  electron  plasma  waves  in  the  plasma  (Tri,relpiece-Gould  modes)  for  the  time 
being,  and  treats  only  the  modified  symmetric  TM  electromagnetic  waves.  The  theory 
predicts  that,  when  driven  by  an  electron  beam,  the  presence  of  a  plasma  in  the  slow  wave 
structure  will  cause  an  increase  in  the  oscillation  frequency,  because  the  dispersion  relation 
for  the  electromagnetic  (TMoi)  wave  is  raised  in  frequency  (see  Fig.  31).  As  can  be  seen 
from  Fig.  32,  the  presence  of  the  plasma  effect  higher  order  modes,  too.  The  temporal 
growth  rates  of  a  high  frequency  mode  approach  those  of  the  fundamental  mode  for  very 
high  plasma  densities  (Fig.  33).  A  simplified  analytical  model  predict  enhancement  in  the 
spatial  growth  rate  and  a  switch  from  absolute  to  convective  instability  at  a  critical  plasma 
density  of  Np  =  2  x  1012  cm-3.  The  predicted  enhancement  in  the  linear  growth  rate  is 
attributed  to  a  decrease  in  the  grope  velocity  of  the  backward  wave  in  the  presence  of  the 
background  plasma. 

Even  though  this  simplified  model  predicts  some  interesting  features  like  enhancement 
of  the  linear  growth  rate,  it  occurs  at  a  plasma  density  which  is  100  times  large  than  those 
used  in  the  experiment.  It  is  clear  that  the  present  linear  theory  will  have  to  be  considerably 
expanded  to  better  explain  the  experimental  results  and  to  be  used  as  a  tool  with  predictive 
capabilities.  The  linear  theory,  of  course,  can  not  be  used  to  analyze  efficiency  enhancement 
which  is  a  non-linear  phenomenon.  The  presence  of  the  plasma  can  cause  a  substantial 
increase  in  the  spatial  growth  rate  of  the  aboslute  instability  (Fig.  34).  For  high  plasma 
density  the  absolute  instability  is  suppressed  and  only  the  convective  instability  remains. 


Wavenumber  kz  (  cm  1) 


Figure  31:  Calculated  dispersion  relation  of  the  plasma  loaded  BWO. 
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Figure  32:  Dispersion  curves  for  three  modes  for  various  plasma  densities  A’p,  =  0, 
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Nb=  S.lExIO’^cm'3  ) 
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Figure  33:  Peak  growth  rates  Im(u)max/2TT  vs  plasma  density. 
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ure  34:  Simplified  theoretical  result  of  the  maximum  spatial  growth  rate  (solid  curve) 
r  !rn i-'l  =  0  normalized  to  the  vacuum  case  vs  plasma  density  ,Vr;  ,\\  —  o.lo  x  IOjO 
’■  The  maximum  temporal  growth  rates  for  Imik)  =  0  normalized  to  the  vacuum  case 

are  shown  by  dashed  curve. 
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2.10  Theory  of  Relativistic  Backward  Wave  Oscillators 

Relativistic  backward  wave  oscillators  (BW'Os)  have  proven  to  be  efficient,  high  power  mi- 

•  'rowave  sources  in  the  centimeter  and  millimeter  wave  range,  capable  of  radiating  hundreds 
■  •I  megawatts  of  power.  One  approach  to  increasing  the  power  generating  capabilities  of 
i  hese  devices  is  to  introduce  plasma  into  the  device.1 

As  a  first  step  in  the  development  of  a  comprehensive  theory  of  plasma  loaded  relativistic 
1.3 WO  operation,  we  have  developed  a  linear  and  nonlinear  theory  for  a  vacuum  relativistic 
BWO  of  a  finite  length.  To  this  end  we  have  derived  a  set  of  reduced  equations  self  consis¬ 
tently  describing  the  slow  evolution  of  the  envelope  of  the  radiation  field  and  the  relativistic 
motion  of  the  particles  under  the  following  assumptions:  (a)  the  electron's  motion  is  one 
dimensional  along  the  z-axis  (axis  of  symmetry)  due  to  a  strong  applied  magnetic  field;  (b) 
t  he  beam  interacts  with  only  one  spatial  harmonic  of  the  electromagnetic  field.  The  model 
takes  into  account  the  finite  reflectivity  of  the  electromagnetic  wave  at  both  boundaries  of 
the  interaction  region  and  allows  large  variation  in  the  longitudinal  velocity  of  the  electron 
beam.  The  model  equations  depend  on  six  normalized  parameters:  the  length.  koL\  the 
initial  beam  velocity,  J.(0)  (corrected  bv  taking  into  account  DC  beam  space  charge  effects); 

•  ne  value  of  the  radiation  group  velocity.  |  J3|;  the  value  of  the  combined  reflection  coefficient 
it  the  boundaries  of  the  RF  structure.  j/?|,  the  current  /  and  the  beam  AC  space  charge 
parameter.  l'sp.  This  is  in  contrast  to  previous  models,2  in  which  radiation  is  assumed  to 
leave  the  interaction  region  without  reflection  and  variations  in  the  longitudinal  velocity  are 
assumed  to  be  small.  Note,  that  in  this  case,  the  number  of  normalized  parameters  can  be 
reduced  to  two:  current  and  A.C.  space  charge  parameter. 

The  normalized  length  is  given  by  k0L  —  Luo/(0*(O)c),  where  c  is  the  speed  of  light,  L 
is  the  length  of  the  structure.  u,-0  is  the  frequency  determined  from  the  intersection  of  the 
dispersion  curve  of  a  corrugated  wave  guide  and  the  beam  dispersion  curve  (see  Fig.  35). 
The  normalized  current  is  given  by  the  following  expression. 


(10) 


where  I  a  =  imc^/q)  is  the  Alfven  current  in  Amperes,  I  is  the  beam  current  in  Amperes,  d 
the  period  of  the  structure,  and  C  is  the  coupling  coefficient. 

The  value  for  C  can  be  obtained  from  the  solution  of  Maxwell’s  equations  in  the  corru¬ 
gated  structure3  without  the  electron  beam  using  the  following  expression 


( ' 


‘2jrdr 

,P 


(|E(r.  r)2|  +  | B( r.  c)|2)  . 


vhere  r,,.|  ~).  (n,)  is  the  wall  radius  (beam  radius)  and  E.  B  and  E,  are  the  electric,  magnetic 
fields  and  longitudinal  component  of  the  electric  field,  respectively.  The  AC  space  charge 
parameters.  J,p  is  the  frequency  of  beam  plasma  oscillations  (including  the  effect  ot  the  beam 
momotryi  normalized  to  the  time  of  flight.  L/n.(0).  The  exact  expression  for  u,',p  will  be 
divert  elsewhere. 

We  calculated  the  start  oscillation  current  for  different  values  of  the  parameters  koL  and 
!!'  and  fixed  the  values  of  the  other  parameters.  The  normalized  start  oscillator  current  is 
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related  to  the  normalized  current  (Eq.  (10))  as  follows 


=  /(M-V  (7o3#(0))  (11) 

Figure  36  displays  the  normalized  start  csciliation  current  as  a  function  of  the  normalized 
length  k0L  for  combined  amplitude  reflection  coefficients  of  0.  0.3  and  0.7.  From  the  plot, 
we  observe  that  for  the  )R\  =  0  case  the  start  current  roughly  has  a  constant  value.  I  ~  7.7. 
For  \R\  ^  0.  on  average,  the  same  scaling  exists  in  addition  to  a  periodic  dependence  on 
kaL.  which  is  due  to  the  interference  between  the  amplified  and  reflected  waves.  The  actual 
value  of  the  start  oscillation  current  can  be  obtained  from  Eqs.  (10)  and  (11)  and  Fig.  36. 
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We  observe  that  by  varying  the  normalized  length  parameter  by  7r/2,  one  can  expect  a 
variation  in  the  start  oscillation  current  of  a  factor  of  three  for  the  R  =  . 7  case.  The 
presence  of  a  low  density  plasma  modifies  the  dispersion  characteristics  of  the  RF  structure. 
Simple  estimates  show  that  a  change  in  the  parameter  k0L  (the  operation  point)  of  order 
t/2  corresponds  to  u>v/uj0  ~  0.3,  where  u ;p  —  (Airnq2 /m)1^  is  the  plasma  frequency  and  n 
is  the  plasma  number  density.  Note,  to  achieve  a  similar  variation  in  the  k0L  parameter 
will  require  about  60%  variation  in  the  voltage  on  the  gun.  The  preceding  results  are  not 
sensitive  to  variations  in  &(0)  and  \0g\.  Figure  37  shows  the  normalized  frequency  shift,  6u, 
at  start  oscillation  versus  k0L  assuming  L  fixed.  Here  6u  =  u(k0L)—uj(k0L  =  30).  Variation 
of  tiui  with  k0L  assuming  the  structure  fixed  occurs  when  beam  energy  is  varied,  and  the 
point  of  intersection  of  the  beam  and  structure  dispersion  curves  illustrated  in  Fig.  35  moves. 
In  the  case  of  \R\  =  0  the  operating  frequency  decreases  with  increasing  k0  reflecting  the 
continuous  changes  in  the  intersection  point.  When  )/?|  ^  0  the  frequency  tends  to  change  in 
discrete  steps.  This  is  a  manifestation  of  the  tendency  of  the  cavity  to  have  eigenfrequencies 
as  the  reflectivity  is  raised.  On  the  average  the  operating  frequency  tracks  the  intersection 
frequency:  however,  with  high  reflectivity  the  frequency  changes  in  steps.  Figure  38  shows 
the  time  history  for  the  electronic  efficiency  for  two  examples:  the  beam  current  is  two 
and  five  times  larger  than  the  start  oscillation  current.  The  normalized  length  k0L  —  30, 
1.(0)  =  .768,  Jg  =  .24,  |/?|  =  0  and  u>sp  =  0.  As  it  can  be  seen,  for  the  larger  current,  the 
radiated  power  has  large  fluctuations  (i.e.,  over-bunch  instability)  and  for  currents  less  than 
'ome  critical  value.  /cr,  the  radiation  power  is  constant.  For  very  large  values  of  k0L  (~  100) 
and  1  R\  =0  the  critical  value  /cr  ~  3/J(art,  which  is  consistent  with  previous  results.2  In  the 
vn.se  of  l/?|  -£  0.  the  Lr  is  a  complicated  function  of  k0L  and  is  shown  on  Fig.  36  together 
with  t he  value  of  /  (Eq.  (10))  at  start  oscillation  for  \flg\  =  0.24  and  d-(0)  =  0.746.  A  similar 
behaviour  of  IrT  versus  k0L  is  anticipated  for  other  values  of  \8g\  and  d.(0).  Preliminary 
estimates  indicate  that  Maryland  BWO  is  operating  in  the  overbunch  regime.  This  part  of 
t  he  work  was  presented  at  a  number  of  conferences  (see  Section  2.12)  and  a  manuscript  is  in 
preparation. 
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Figure  35:  Dispersion  diagram  of  the  BWO. 


igure  36:  Normalized  start  oscillation  current.  /,  vs.  normalized  length.  knL,  for  a  BWO 
with  R= 0.  .3.  and  .7  solid  line,  dashed  line,  and  large  dashed  line,  respectively. 
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Figure  37:  Normalized  frequency  shift.  6u,  versus  normalized  length  koL  (fixed  L)  at  the 
start  oscillation  current  for  BWO  with  R  =  0,  .3  and  .7  solid  line,  dashed  line,  and  large 

dashed  line,  respectively. 


Normalized  Time  (T=tv,/L  ) 


Figure  38:  Efficiency  of  relativistic  BWO  versus  normalized  time  with  R  =  0  for 
/  =  o  l, tan  l solid  line)  and  1  =  2  l,tart  dashed  line. 
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2.11  Open  Research  Topics 

We  achieved  considerable  progress  on  both  the  experimental  and  theoretical  aspects  of  high 
power  microwave  generation  using  plasma  loaded  BWO.  The  work  is  not  complete,  though, 
and  additional  research  is  needed  in  order  to  take  full  benefit  of  fhe  innovative  approach  of 
plasma  microwave  electronics.  Experimentally,  the  following  subjects  need  to  be  addressed 
i  partial  list ): 

(a)  operation  at  currents  well  above  the  vacuum  space  charge  limiting  value 

(b)  measurement  of  the  starting  currents  in  the  presence  of  plasma 

(c)  optimization  with  respect  to  interaction,  efficiency  and  pulse  length 

(d)  improved  diagnostics  to  probe  the  plasma  processes,  and  improved  plasma  gener¬ 
ation  techniques 

(e)  investigation  of  overmoded  systems 

Theoretically,  our  work  needs  to  be  expanded  in  order  to  gain  a  better  understanding  of 
the  three-component  system  ( beam/ plasma/ wave),  finite  magnetic  field  and  space  charge 
effects,  and  interaction  in  overmoded-multiwave  systems.  It  is  expected  that  the  combined 
experimental  and  theoretical  efforts  will  lead  to  major  achievement  in  this  area  and  in  plasma 
microwave  electronics. 
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Absolute  Instability  for  Enhanced  Radiation  from  a  High-Power  Plasma-Filled 

Backward-Wave  Oscillator 
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The  linear  theory  of  electromagnetic  radiation  from  a  backward-wave  oscillator  with  a  plasma-filled, 
sinusoidally  corrugated  waveguide  driven  by  a  relativistic  electron  beam  has  been  derived  and  analyzed 
numerically.  The  presence  of  plasma  can  cause  a  substantial  increase  in  the  spatial  growth  rate  of  the 
absolute  instability.  For  high  plasma  densities,  however,  the  absolute  instability  is  suppressed  and  only 
the  convective  instability  remains.  The  predicted  radiation  enhancement  can  be  attributed  to  a  decrease 
in  the  group  velocity  of  the  backward  wave  in  the  presence  of  a  plasma. 

\ 

PACS  numbers;  S2.40.-w,  4l.80.Ee,  52.2S.Sw,  85.10.Ji 


In  recent  years,  a  dramatic  resurgence  in  research  ac¬ 
tivities  on  relativistic  electronics  for  high-power  mi¬ 
crowave  sources  has  been  taking  place  all  over  the  world, 
driven  by  the  development  of  new  technologies  and  the 
requirements  of  present  and  future  applications. 1  Gyro- 
trons  and  free-elcctron  lasers  are  examples  of  such  new 
devices.  The  backward-wave  oscillator  (BWO)  driven 
by  an  intense  relativistic  electron  beam  is  another  mi¬ 
crowave  source  on  which  research  has  been  conducted 
primarily  in  the  U.S.  and  U.S.S.R.  for  the  last  fifteen 
years.  :*'°  Previously,  gas-tilled  BWOs  were  observed  to 
yield  significant  microwave  output  power  at  A-band  fre¬ 
quencies.3'9  Recently,  however,  considerable  increases  in 
both  the  output  power  and  the  efficiency  of  a  BWO  due 
to  plasma  injection  from  an  external  plasma  gun  have 
been  observed.10  Lin  and  Chen"  have  recently  reported 
numerical  simulations  of  the  experiments,  although  they 
assumed  an  artificial  periodic  boundary  condition  in  the 
axial  direction  which  did  not  correspond  to  the  experi¬ 
ments.  They  reported  that  a  parametric  beam-plasma 
wave  interaction  could  not  be  the  dominant  mechanism 
for  the  efficiency  enhancement  and  that  the  enhancement 
could  be  attributed  to  a  decrease  in  the  phase  velocity  of 
the  most  unstable  beam  mode  in  the  presence  of  a  back¬ 
ground  plasma.  In  order  to  understand  analytically  the 
physical  reason  for  the  observed  enhancement,  we  have 
developed  a  linear  theory  of  the  beam-plasma  interaction 
and  have  analyzed  numerically  the  absolute  instability12 
in  plasma-filled  BWOs. 

First,  we  consider  an  infinitely  long  axisymmetric 
waveguide,  i.e.,  slow-wave  structure,  whose  wall  radius 
varies  with  the  axial  coordinate  z  according  to  R(z) 
“■/?o+hcos(koz),  where  ko"2ir/zo.  Here  Ro,  h,  and  zo 
are,  respectively,  the  average  waveguide  radius,  the  am¬ 
plitude  of  corrugation,  and  its  period.  A  uniform,  cold, 


and  collisionless  plasma  with  density  Np  is  assumed  to  be 
present,  and  a  beam  with  uniform  electron  density  /V*,  a 
longitudinal  velocity  u,  and  radius  Rk<R0-h  is  present 
in  the  slow-wave  structure.  An  infinitely  strong  external 
magnetic  field  which  confines  electron  motion  to  be 
strictly  along  the  field  lines  is  applied  in  the  axial  direc¬ 
tion  z.  We  have  obtained  the  linear  dispersion  relation 
Dlk,ta)  “0  for  the  symmetric  TM  modes  with  field  com¬ 
ponents  £,.  E„  and  if#/3  Here,  D  is  the  value  of  the 
determinant  of  a  matrix  for  which  the  element  of  the 
mth  column  and  nth  row,  D^,  is  given  by 

Dmm-[\  +(n-m)QMKHCL+LnCD ,  (1) 

""  ’  ?-o  2*+l"-m|9!(q  +  |n-m|)!  ’ 

K,  -  YnJ^HS)N\{YH8) -X„JxU(n8)N^Y,8) , 


L,  -X'Jo(Y'8)JM)  -  YJoiX'SMiiYnS) . 


Pip _ Oil 

a>2  y3(o>  — /t„i>)2 


(2) 


Qn~ 


kok„ 


w2/c2  —  ki 


R  o 


A  similar  expression  for  Cm,  holds,  with  J o  replaced 
by  No  in  CJmH  given  by  (2).  Here,  we  have  assumed  a 
phase  factor  expl  — i(a»  —  k„r)]  for  time  and  spatial 
variations  of  the  rf  electric  fields,  and  the  other  notations 
are  standard.  We  confine  our  analysis  here  to  the  case  of 
the  high-frequency  electromagnetic  TMoi  mode  with 
at  >  Up.  Our  dispersion  relation  ( 1 )  tends  to  previously 
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obtained  results'*'3  in  the  limit  of  small  a. 

The  following  parameters9,10  are  adopted  unless 
specific^  otherwise;  the  relativistic  factor  y  —  2.23,  Ro 
“1.445  cm,  R),  “0.9  cm,  h  ”0.445  cm,  co“l  -67  cm, 
and  /V„ -2.09x10"  cm-3.  The  plasma  density  Np  is 
considered  to  be  a  variable.  For  these  parameters,  we 
must  calculate  a  9x9  determinant  equation,  with  the 
•~rms  up  to  the  order  of  a10  in  (2),  to  obtain  1%  numeri¬ 
cal  accuracy  in  periodicity  for  two  periods  of  wave  num¬ 
ber  ko . 13 

The  spatial  growth  rate  is  a  more  practical  measure  of 
the  strength  of  instabilities  than  is  the  temporal  growth 
rate,12  as  will  be  discussed  later  in  detail.  There  can  be 
many  roots  of  D(k,u>)  “0  that  have  complex  k  — k,+ikj 
for  a  real  co.  In  order  to  sort  out  the  spatially  growing 
waves  of  convective  and  absolute  instabilities  from  other 
roots  of  evanescent  or  stable  waves,  we  deform  the  path 
of  integration  in  the  inverse  Laplace  transformation  in 
the  complex  <u  “  <u, -Htu,  plane  downward  to  the  real 
axis.  If  k,  has  different  signs  when  a,  takes  a  large  posi¬ 
tive  value  and  zero,  then  the  wave  is  convectively  unsta¬ 
ble:  otherwise,  it  is  an  evanescent  wave.  If  the  deforma¬ 
tion  of  the  path  of  integration  towards  the  real  axis  is 
prevented  by  merging  of  two  roots  (saddle  point)  for 
some  positive  to,  from  both  sides  of  the  k  half  planes, 
then  the  instability  is  absolute  which  corresponds  to  a 
BWO.  In  general,  to  and  k  at  the  saddle  point  are  both 
complex1214  (*,  +  ui,  “ai„ +iw„),  and  the  ab¬ 

solute  instability  causes  electromagnetic  radiation  from 
an  infinitely  long  slow-wave  structure  even  if  a  feedback 
mechanism  does  not  exist.  The  spatial  growth  rate  k„ 
and  the  real  part  k„  of  the  wave  number  at  the  saddle 
point  versus  A'„  are  plotted  in  Fig.  1  (a),  while  the  tem¬ 
poral  growth  rate  w„/2/r  and  the  oscillation  frequency 
to,,/ 2; r  are  shown  in  Fig.  1  (b).  As  /V,,  increases  the  sad¬ 
dle  point  goes  down,  and  the  spatial  growth  rate  is 
enhanced.  As  the  point  moves  into  the  to,  <  0  region  as 
shown  in  Fig.  Kb),  the  absolute  instability  disappears 
suddenly  and  only  the  convective  instabilities  continue  to 
exist  for  exceedingly  large  Np.  These  convective  insta¬ 
bilities  [corresponding  to  a  plasma  traveling-wave  tube 
(TWT)l  cannot  give  rise  to  any  radiation  unless  a  feed¬ 
back  mechanism  exists  in  the  slow-wave  structure.  The 
maximum  spatial  growth  rates  k,m  for  ai,  “0  are  a  mea¬ 
sure  of  the  amplification  in  the  plasma  BWO  and  are 
shown  in  Fig.  1(a).  The  group  velocity,15  approximately 
given  by  co„lk„,  and  the  phase  velocity  <o„/k„  vs  Np  are 
shown  in  Fig.  2  for  an  infinitely  long  (£.“<»)  plasma- 
loaded,  slow-wave  structure.  The  absolute  value  of  the 
group  velocity  decreases  with  Np  while  the  phase  velocity 
slightly  increases  with  Np.  The  group  velocity,  however, 
changes  direction  and  becomes  positive  for  plasma  densi¬ 
ty  Np  >  1 .6  x  1 0 12  cm"3  as  shown  in  Fig.  2  (solid  cir¬ 
cles).  in  which  case  the  system  works  as  a  plasma  TWT 
(a)i  <  0),  because  the  wave  is  still  growing  spatially. 

In  our  linear  theory,  the  output  radiation  is  propor¬ 
tional  to  exp(2ci>„<).  where  /  is  the  interaction  time  be- 


PLASMA  DENSITY  Np(xl010anJ  ) 


FIG.  I.  (a)  The  spatial  growth  rale  k„,  the  real  part  k,„  at 
the  saddle  jxtint.  and  the  maximum  spatial  growth  rate  for 
lm(w)”0,  k,m,  vs  plasma  density  Nn.  (b)  Temporal  erowth 
rate  w„/2jr,  oscillation  frequency  <o„/2/c.  at  the  saddle  point, 
and  plasma  frequency  top/2it  vs  Np. 


tween  the  beam  and  the  backward  wave.  In  real  experi¬ 
ments,  the  interaction  time  f»l  nsec  is  given  by  the 
ratio  of  a  finite  L  and  the  group  velocity.  Then. 
exp(2«/,f)  “exp(2£„L)»  1.  The  increase  in  k„  with 
N,,  as  shown  in  Fig.  1  (a)  may  thus  enhance  the  radia¬ 
tion.  This  is  the  reason  why  the  spatial  growth  rate  k„  is 
a  more  practical  and  critical  measure  of  the  instabilities 
than  the  temporal  growth  rate  a>„/2/r  for  a  slow-wave 
structure  with  finite  axial  length.  The  predicted  en¬ 
hancement  can  be  attributed  to  an  increase  in  the  in¬ 
teraction  ume,  i.e„  a  decrease  in  the  group  velocity  of 
the  backward  wave.  As  shown  in  Fig.  1(a),  our  result 
predicts  linear  enhancement  for  values  of  Np  1  order 
greater  than  that  observed  in  Ref.  1 1. 

In  order  to  elucidate  the  physical  meaning  of  the  radi¬ 
ation  enhancement,  the  maximum  spatial  growth  rate 
k,m  of  the  plasma  BWO  in  the  limit  of  small  a  has  been 
analyzed  in  a  qualitative  way.  The  procedure  used  here 
is  not  very  different  from  that  in  Ref.  1 1.  A  2x2  deter¬ 
minant  equation  D(k.co)  “0,  with  —  l<m,n<0,  is 
used,  and  only  the  terms  up  to  order  a  are  taken  into  ac¬ 
count  in  (2).  We  assume  that  the  beam  term  in  A,  [see 
Eq.  (2)1  is  small,  and  the  spatial  growth  rate  (o>,  “0)  at 
the  intersection  of  the  beam  space-charge  wave  co~kv 
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FIG.  2.  Phase  and  group  velocities  of  the  plasma-loaded 
slow-wave  structure  in  the  case  of  L  “»  vs  N shown,  respec¬ 
tively.  by  solid  and  open  circles.  Those  for  a  slow-wave  struc¬ 
ture  with  finite  axial  length  Z.“I6.7  cm  are  shown  by  thin 
lines. 


and  the  backward  wave  Y\  “2.405  is  calculated.  The 
final  dispersion  relation  results  in  a  quartic  equation,  it 
must  be  emphasized  here  that  the  previous  analytical  re¬ 
sult5  for  y»  I  gives  no  enhanced  growth  rates  for  the 
parameters  in  the  experiments910  with  y“2.23.  The 
spatial  growth  rate  versus  calculated  from  the  quartic 
equation  is  shown  in  Fig.  3.  A  resonant  increase  in  the 
spatial  growth  rates  is  found  and  the  result  is  qualitative¬ 
ly  similar  to  k,m  in  Fig.  1(a).  If  one  ignores  the  fourth- 
order  term  in  the  quartic  equation  assuming  that  it  is 
small,  the  peak  in  Fig.  3  goes  to  infinity  at  kr  “fco  and  at 
the  cutoff  frequency,  where  the  group  velocity  of  the 
backward  wave  is  zero.  The  decrease  in  the  group  veloc¬ 
ity  enhances  the  interaction  between  the  backward  wave 
and  the  beam.  This  is  the  physical  reason  for  the 
enhanced  radiation  from  the  plasma  BWO  predicted  in 
the  present  analysis.  The  temporal  growth  rate  ()c,  “0) 
which  corresponds  to  <u/,/2jt  in  Fig.  1  (b)  has  also  been 
calculated  in  the  same  2x2  determinant  equation,  and 
the  result  is  shown  by  a  dashed  line  in  Fig.  3.  It  should 
be  noted  that  the  temporal  growth  rate  is  zero  for 
Np  >  l.7x  1012  cm  ~3. 

So  far  we  have  treated  an  infinitely  long  slow-wave 
structure.  For  a  finite-length  device,  different  modes 
may  be  linked.  In  our  analysis  we  considered  four 
modes — two  structure  modes  (forward  and  backward 
TMoi)  and  two  beam  space-charge  modes  (four  roots  in 
the  dispersion  relation).  The  oscillation  is  affected  by 
the  reflection  coefficients  R\  and  Ri  at  both  ends. 14  The 
oscillation  does  not  occur  at  the  k,  of  the  saddle  point, 
but  at  a  growing  root  k+  near  k,,  where  k+  is  found 


PLASMA  DENSITY  Np(xl010cm'3  ) 


FIG.  3.  Simplified  theoretical  result  of  the  maximum  spatial 
growth  rate  (solid  curve)  for  Im(o>)  “0  normalized  to  the  vac¬ 
uum  case  vs  plasma  density  N N*“5.15x  IO'#cm“J,  a  “0.1 
and  other  parameters  are  the  same  as  those  in  Figs.  1  and  2. 
The  maximum  temporal  growth  rates  for  Im(fc)  “0  normalized 
to  the  vacuum  case  are  also  shown  by  the  dashed  curve. 

from  the  following  simultaneous  equations: 

D(JU,a>)“0,  D(k  +  +  f,o>)  “0  , 

e-k--k+--2xN/L~i\n(R)/L  . 

and  A  — ±  1.  ±2 . 

Here,  k  -  is  the  evanescent  root  propagating  in  the  oppo¬ 
site  direction  to  k  ♦  in  order  to  feed  back  the  oscillation 
energy  to  the  original  position.  In  this  work,  we  neglect 
mode  conversion  at  both  ends  and  consider  only  the 
growing  wave  with  highest  growth  rate  k  +  and  the 
evanescent  wave  with  smallest  damping  rate  k  -  because 
the  main  contribution  comes  from  these  two  waves. 
Note  that,  to  have  oscillation,  in  addition  to  o>,  >  0,  the 
spatial  growth  rate  Im(/fc  +  )  must  surpass  the  spatial 
damping  rate  Im(jfc-),  because  £<  1.0.  This  imposes 
an  additional  restriction  to  backward-wave  oscillation. 
We  do  not  find  k+  in  the  case  of  convective  instabilities. 
The  phase  and  group  velocities  for  Z-  “  1 6.7  cm  and 
R  “1.0  are  shown  in  Fig.  2.  The  thin  line  for  phase  ve¬ 
locity  denotes  the  physical  solutions  which  change  from 
N  “  —  I  to  A  “  I  as  Np  increases.  The  plasma  BWOs 
will  not  operate  (o>i  <  0)  for  L  and  R  less  than  the  criti¬ 
cal  values13  for  given  Nr  and  Np. 

In  conclusion,  we  present  a  small-signal-gain  calcula¬ 
tion  for  a  plasma-loaded  slow-wave  structure  immersed 
in  an  infinitely  large  axial  magnetic  field  for  frequencies 
at  >  atp.  It  was  found  that  the  presence  of  the  plasma 
modifies  both  the  dispersion  relation  and  the  nature  of 
the  vacuum  linear  instability.  For  low  plasma  densities, 
the  instability  is  absolute.  For  large  plasma  densities, 
the  group  velocity  changes  sign,  in  which  case  the  insta¬ 
bility  became  convective.  The  present  analysis  predicts  a 
considerable  increase  both  in  the  spatial  growth  rate  and 
in  the  oscillation  frequency  with  increased  Np  (Fig. 
1(b)).  The  saturated  levels  and  efficiencies11  of  plasma 
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BWOs  are  beyond  the  scope  of  the  present  linear 
analysis.  Nonlinear  calculations  are  underway. 
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Abstract— ka  experimental  demonstration  of  a  strong  enhancement 
of  the  interaction  efficiency  in  a  high-power  relativistic  backward-wave 
oscillator  when  plasma  is  injected  is  presented.  Controlled  plasma  in¬ 
jection  enhances  the  interaction  efficiency  over  the  vacuum  case  by  a 
factor  of  up  to  8  to  a  value  of  about  40%.  A  linear  theory  of  electro¬ 
magnetic  wave  generation  in  plasma-loaded  corrugated  wall  resonators 
is  reviewed.  A  number  of  physical  mechanisms  are  considered  to  ac¬ 
count  for  the  enhanced  interaction,  including  two  variations  of  a  three- 
wave  interaction  involving  the  electron-beam  slow  space-charge  wave, 
the  slow  electromagnetic  waves  in  the  structure,  and  the  quasi-electro¬ 
static  waves  in  the  plasma. 

I.  Introduction 

HE  FIELD  of  plasma  electronics  dates  back  to  1949 
when  several  authors  [1]  discussed  the  excitation  of 
electromagnetic  waves  in  a  plasma  by  an  electron  beam 
by  virtue  of  stimulated  Cerenkov  radiation.  Since  1949 
this  field  of  research  has  expanded  to  include  beam  heat¬ 
ing  of  plasmas,  collective  acceleration  of  charged  parti¬ 
cles  by  plasma  fields,  and  the  excitation  and  amplification 
of  electromagnetic  waves.  This  last  area  has  been  referred 
to  in  the  literature  as  “plasma  microwave  electronics” 
[21. 

The  introduction  of  plasma  into  vacuum  microwave  de¬ 
vices  can  have  several  beneficial  effects.  It  is  expected 
that  plasma  microwave  devices  could  be  tuned  in  fre¬ 
quency  by  exercising  control  over  the  plasma  density  up 
to  the  millimeter  wave  regime.  Finally,  extremely  high- 
power  microwave  devices  may  be  realized  by  using  beam 
currents  well  above  the  space-charge-limiting  current  in 
vacuum.  These  currents  may  be  realized  by  using  a  plasma 
to  neutralize  the  beam  space  charge. 

The  practical  realization  of  plasma  microwave  devices 
involves  serious  difficulties  which  have  yet  to  be  fully 
overcome.  These  difficulties  include  the  generation  of 
uniform,  steady-state,  highly  ionized  plasmas,  high  noise 
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levels  which  may  have  an  adverse  effect  on  the  operation 
of  microwave  amplifiers  depending  on  the  application  en¬ 
visioned,  the  lack  of  a  fully  developed  linear  and  nonlin¬ 
ear  theory  of  plasma  devices,  and  the  problem  of  coupling 
the  electromagnetic  radiation  both  into  and  out  of  the 
plasma. 

In  particular,  the  problem  of  coupling  the  electromag¬ 
netic  radiation  to  plasma  waves  has  proved  to  be  quite 
complicated  and  must  be  resolved  before  truly  effective 
plasma  microwave  devices  may  be  realized.  The  difficulty 
lies  with  the  slow  longitudinal  oscillations  which  are  ex¬ 
cited  in  the  plasma.  For  nonrelativistic  beams,  these  os¬ 
cillations  have  very  slow  phase  velocities  (fp^  «  c) 
and  behave  like  quasi-electrostatic  oscillations.  As  a  con¬ 
sequence,  they  are  mostly  trapped  in  the  plasma.  Recent 
efforts  have  sought  to  reduce  the  losses  during  the  extrac¬ 
tion  of  energy  from  a  plasma  by  using  relativistic  electron 
beams.  Under  these  conditions  the  wave  excited  in  the 
plasma  may  have  a  phase  velocity  vphtx  ~  c,  allowing  it 
to  be  more  easily  extracted  from  the  plasma. 

The  purpose  of  this  work  is  to  study  the  feasibility  of 
generating  very  high-power  microwave  radiation  by  using 
relativistic  electron  beams  in  plasma-loaded  microwave 
devices.  Over  the  years  a  wide  variety  of  high-power  mi¬ 
crowave  (HPM)  vacuum  devices  have  been  studied,  in¬ 
cluding  the  well-known  magnetron,  klystron,  and  back¬ 
ward  wave  oscillator,  as  well  as  newer  devices  such  as 
the  gyrotron,  the  free-electron  laser,  and  the  virtual  cath¬ 
ode  oscillator  [3],  All  of  these  sources  have  one  thing  in 
common— they  are  driven  by  an  intense  unneutralized 
electron  beam  which  interacts  unstably  in  high  vacuum 
with  an  electromagnetic  wave,  leading  to  the  conversion 
of  the  kinetic  energy  of  the  beam  into  electromagnetic  ra¬ 
diation.  The  power  levels  available  from  such  devices 
have  grown  by  an  order  of  magnitude  every  decade  since 
1940.  For  example,  advanced,  large  diameter,  over- 
moded  backward  wave  oscillators  (like  the  multiwave 
Cerenkov  and  multiwave  diffraction  generator)  have'  re¬ 
cently  achieved  efficient  operation  at  power  levels  of  about 
15  GW  [4]  at  a  wavelength  of  3  cm  and  4.5  GW  [5]  at  a 
wavelength  of  6.5  mm,  and  operation  at  a  wavelength  of 
2  mm  has  been  demonstrated  [6],  The  thrust  of  this  work 
is  to  evaluate  efficiency  enhancement  in  a  small  diameter, 
relativistic  backward-wave  oscillator  using  a  different  ap¬ 
proach— that  of  plasma  loading. 
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Backward-wave  oscillators  (BWO's)  are  ideal  candi¬ 
dates  for  the  evaluation  of  plasma  effects  on  microwave 
generation  because  they  are  relatively  simple  devices, 
providing  fairly  effective  conversion  of  electron-beam  en¬ 
ergy  into  radiation,  and  can  be  easily  filled  with  plasma. 
Initial  studies  [7],  [8]  have  utilized  background  plasmas 
produced  by  the  electron-beam  impact  ionization  of  a  low- 
pressure  neutral  gas  background  and  have  demonstrated 
enhanced  microwave  output  power.  In  recent  experiments 
high-power  microwave  radiation  was  observed  in  a  rela¬ 
tivistic  backward-wave  oscillator  which  was  externally 
filled  with  a  highly  ionized  plasma  {9].  Controlled  plasma 
injection  enhanced  the  interaction  efficiency  relative  to  the 
vacuum  case  by  as  much  as  a  factor  of  8,  resulting  in  an 
efficiency  of  about  40%. 

Section  II  of  this  paper  reviews  wave  propagation  in  a 
plasma-filled  smooth-walled  waveguide;  Section  III  con¬ 
tains  a  discussion  of  the  space-charge-limiting  current  in 
plasma-loaded  waveguides.  A  linear  theory  of  plasma- 
filled  corrugated  structures  is  reviewed  in  Section  IV. 
Section  V  summarizes  the  results  from  an  experimental, 
plasma-loaded,  high-power  BWO  that  demonstrated  a 
strong  enhancement  in  the  interaction  efficiency,  and  Sec¬ 
tion  VI  contains  a  discussion  of  the  results,  including  such 
questions  as  possible  coupling  mechanisms,  mode  stabil¬ 
ity,  and  future  trends  in  plasma  microwave  electronics. 

II.  Review  of  Wave  Propagation  in  Plasma-Filled, 
Smooth-Walled,  Cylindrical  Waveguides 

Plasma  effects  in  conventional  microwave-generating 
devices  can  usually  be  neglected  since  the  plasma  fre¬ 
quency  of  the  background  gas  is  much  smaller  than  the 
plasma  frequency  of  the  electron  beam.  Typically,  wp/u>b 
<  10~2,  where  up  =  (eNp/m,t0)'/2  is  the  background 
plasma  frequency,  and  w*  =  (eNh/mrt0)>/2  is  the  beam- 
plasma  frequency.  Recent  theoretical  studies,  however, 
have  predicted  that  the  presence  of  a  plasma  in  high-power 
microwave  devices  may  lead  to  enhanced  performance, 
attracting  renewed  scientific  interest  [2],  [  10]— [  1 6). 

One  of  the  characteristics  of  a  fully  ionized  plasma  is 
its  ability  to  support  electric  fields  of  almost  arbitrarily 
targe  amplitude.  A  small  local  deviation  from  neutrality 
in  an  otherwise  neutral  plasma  can  give  rise  to  very  large 
local  fields.  As  an  example,  consider  a  uniform  plasma 
with  an  average  density  of  1014  cm'3.  A  1%  fluctuation 
in  the  local  electron  density  over  a  distance  of  1  mm  pro¬ 
duces  a  local  electric  field  of  -  106  V/cm.  This  capabil¬ 
ity  of  sustaining  very  high  electric  fields  makes  plasma 
attractive  for  use  in  high-power  density  microwave  gen¬ 
erators. 

The  presence  of  a  plasma  in  a  waveguide  changes  the 
characteristics  of  the  guide  and  adds  new  slow  modes  of 
propagation.  For  the  purpose  of  this  discussion,  the  term 
“plasma"  will  be  used  to  describe  a  fully  ionized  gas 
which,  in  the  absence  of  external  disturbances,  is  electri¬ 
cally  neutral.  It  is  further  assumed  that  the  ions  in  the 
plasma  are  stationary,  and  that  the  electrons  have  no  ther¬ 
mal  or  random  velocities  and  undergo  no  collisions.  In 


other  words,  we  will  consider  an  ideal  electron  plasma. 
The  term  “plasmaguide  mode"  will  be  used  to  denot  ny 
of  the  additional  waveguide  modes  that  exist  due  to  the 
presence  of  the  plasma  inside  of  the  guide.  These  are  the 
slow-wave  modes  on  a  nondrifting,  ion-neutralized  plasma 
column  [17J. 

In  order  to  generate  microwave  radiation,  we  are  inter¬ 
ested  in  waves  that  can  cause  bunching  in  an  axially 
streaming  relativistic  electron  beam.  Therefore  only 
waves  with  axial  components  E,  will  be  considered.  There 
are  two  such  families  of  waves;  The  first  is  the  familiar 
TM/,  family  of  electromagnetic  modes  which  will  be 
somewhat  modified  by  the  presence  of  the  plasma,  and 
the  second  is  a  family  of  plasmaguide  inodes 

As  was  demonstrated  by  Trivelpiece  and  Gould  (T-G) 
[17J,  the  dispersion  relation  for  a  plasma-loaded,  smooth- 
walled,  cylindrical  waveguide  immersed  in  an  infinitely 
large  axial  magnetic  field  is  given  by: 


CO 


^  =  k2.  + 


(pjRoV 

i  2,2 

1  -  Wp/uJ 
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where  u  is  the  angular  frequency;  k.  is  the  axial  wave 
number;  c  is  the  speed  of  light  in  vacuum;  R0  is  the  wave¬ 
guide  radius;  and  p,,  is  the  rth  root  of  the  /th  order  Bessel 
function  of  the  first  kind. 

If  the  operating  frequency  of  the  device  is  greater  than 
the  plasma  frequency  (w  >  u>p),  the  electromagnetic 
waves  in  the  plasma-loaded  waveguide  are  very  similar  to 
the  TM  waves  in  a  smooth-walled  empty  waveguide,  as 
can  be  seen  from  (1).  The  only  difference  is  that  the  cut¬ 
off  frequency  for  the  plasma-filled  guide  is  increased  rel¬ 
ative  to  the  empty  waveguide  by: 


u 


2 

CO 


(2) 


This  shift  in  the  cutoff  frequency  can  be  seen  in  Fig.  1(b), 
where  the  dispersion  curve  for  the  lowest-order  symmet¬ 
ric  electromagnetic  mode  (labeled  “TM0i")  is  displayed. 
In  addition,  the  presence  of  the  plasma  in  the  waveguide 
allows  k2  to  take  on  positive  real  values  for  <o  <  cop,  thus 
giving  rise  to  a  propagating  wave.  Fig.  1(b)  also  displays 
plots  of  several  plasma-guide  (T-G)  modes.  The  geome¬ 
try  of  the  oblem  is  shown  in  Fig.  1(a). 

Thus  the  lower  branches  in  Fig.  1(b)  represent  a  family 
of  plasma-guide  modes  having  a  high-frequency  cutoff  at 
up  and  no  low-frequency  cutoff.  One  of  the  interesting 
features  c  he  plasma-guide  modes  is  that  the  high-fre¬ 
quency  cutoff  is  independent  of  the  geometry  and  depends 
only  on  the  plasma  frequency.  This  is  in  contrast  with  the 
electromagnetic  waveguide  modes  whose  cutoff  frequen¬ 
cies  depend  intimately  on  the  geometry.  Another  differ¬ 
ence  between  the  plasma-guide  and  electromagnetic 
modes  is  the  fact  that  for  operating  frequencies  within  the 
plasma-guide  mode  band  (0  <  w  <  up),  all  of  the  plasma¬ 
guide  modes  ( /  2:  0,  v  s  1 )  will  propagate  simulta¬ 
neously  if  excited.  In  contrast,  at  frequencies  above  the 
plasma  frequency  the  number  of  electromagnetic  wave- 
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Fig.  I.  (a)  A  smooth-wall  plasma-loaded  waveguide  immersed  in  an  infi¬ 
nitely  large  guiding  magnetic  field  (model  of  Section  11).  (b)  The  dis¬ 
persion  relation  for  a  plasma-loaded  smooth-wall  waveguide  having  a 
radius  of  R„  =  1.44  cm.  A  single  TMoi  mode  wun  and  without  a  plasma, 
as  well  as  several  plasma-guide  (T-G)  modes,  are  shown. 


guide  modes  that  can  propagate  at  any  given  frequency  is 
finite,  though  increasing  with  frequency.  The  principal 
feature  of  plasma-guide  propagation  is  that  a  plasma  col¬ 
umn  can  support  modes  of  wave  propagation  below  the 
plasma  frequency,  even  in  the  absence  of  electron  drift 
motion. 

As  was  pointed  out  earlier,  the  plasma-guide  modes  are 
electromechanical  in  nature.  The  wave  propagation  re¬ 
sults  from  the  interchange  between  the  kinetic  energy  of 
the  electrons  and  the  stored  energy  in  the  electric  field.  In 
contrast,  for  the  case  of  electromagnetic  waves  in  an 
evacuated  waveguide,  the  wave  propagation  results  from 
the  interchange  of  the  electric  and  magnetic  stored  en¬ 
ergy.  Nevertheless,  there  are  some  similarities  between 
the  two  families  of  waves  in  terms  of  the  field  compo¬ 
nents,  mode  structure,  and  power-carrying  capabilities. 
Table  I  compares  the  field  components  and  power  carried 
by  the  two  classes  of  waves  for  the  case  of  azimuthal  sym¬ 
metry  (/  =  0). 

The  two  families  of  modes  propagate  in  two  different 
frequency  ranges.  While  the  T-G  (plasma-guide)  modes 
propogate  below  the  plasma  frequency,  the  electromag¬ 
netic  modes  can  only  propagate  at  frequencies  above: 


w  > 


(3) 


Using  Table  I,  a  number  of  informative  plots  for  the 
T-G  modes  may  be  generated.  For  example,  Fig.  2  shows 
the  electric  field  configuration  for  the  two  lowest-order 
T-G  modes  (T-G(0.  n  and  T-G(0,2))-  The  power  carried 
by  these  same  modes  as  a  function  of  frequency  is  plotted 
in  Fig.  3,  where  a  plasma  frequency  of  4  GHz  and  a  max¬ 
imum  electric  field  amplitude  of  100  kV /cm  has  been  as¬ 
sumed.  From  this  figure  it  is  clear  that  the  presence  of  the 
plasma  inside  the  waveguide  does  not  lower  the  power¬ 
handling  capabilities  of  the  system.  As  the  wave  fre¬ 


TABLE I 

Field  Components  and  Axial  Power  Flow  for  Azimuthalli 
Symmetric  Electromagnetic  (TM0, )  and  Plasma-Guide  (T-G)  Modes 
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k,  =  pujRa.  Cut-off  wavenumber  for  the  specific  mode  of  interest,  a 
=  £. (r  =  0):  Electric  field  amplitudes  on  axis. 
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Fig.  2.  Electric  field  configuration  for  the  two  lowest-order.  symmetric- 
plasmaguide  modes,  (a)  T-G,0  (b)T-G,0.;,.  Both  were  calculated  for 

u>p  =  2r  ■  4  GHz,  R0  =  1.44  cm  at  a  frequency  of  2  GHz. 

quency  of  the  T-G  mode  approaches  the  plasma  fre¬ 
quency,  the  wave  becomes  increasingly  electrostatic  in 
nature,  its  power-carrying  capability  diminishes,  and  its 
group  velocity  du/dk  decreases  to  zero.  Fig.  4  shows  the 
radial  distribution  of  the  three  field  components  £.,  Er, 
and  H9  for  the  lowest  order  ( /  =  0,  v  -  1 )  symmetric 
plasma-guide  mode  for  various  frequencies.  As  before,  as 
u  approaches  u>p  =  4  GHz  (in  this  example),  both  H$  and 
E,  decrease  to  zero.  Only  the  electrostatic  component  E. 
remains  unchanged. 

III.  Space-Charge-Limiting  Current  in  Plasma- 
Loaded  Waveguides 

In  essence,  a  backward-wave  oscillator  consists  of  an 
electron  beam  confined  radially  by  a  strong  longitudinal 
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Fig.  3.  Calculated  power-carrying  capabilities  of  the  two  lowest-order 
symmetric  plasma-guide  modes  versus  frequency  (u,  =  2i  •  4  GHz). 


Fig.  4.  Calculated  radial  dependence  of  £.,  E„  and  H,  for  the  lowest-or¬ 
der,  symmetric  plasma-guide  mode  at  few  frequencies  (wp  =  2r  ■  4 
GHz,  “  1.44  cm). 


I  z0*  2*7k0i 

Fig.  5.  The  slow  wave  structure  and  beam  model  of  Section  ill. 

the  electromagnetic  field  of  the  struct  -.  The  backward- 
wave  oscillator  is  so  named  because,  at  the  point  of  beam 
and  electromagnetic  structure-mode  resonance,  the  struc¬ 
ture  mode  has  a  negative  group  velocity,  resulting  in  the 
propagation  of  electromagne-  c  energy  backward  along  the 
beam— i.e.,  the  Poynting  vector  points  antiparallel  to  the 
beam  velocity. 

As  indicated  in  Section  I,  recent  theoretical  studies  have 
predicted  that  the  introduction  of  a  plasma  into  vacuum 
HPM  devices  may  lead  to  greatly  enhanced  performance 
[2],  [  10]— [  16] .  For  the  purposes  of  the  discussion  of  the 
space-charge-limiting  current,  we  will  approximate  the 
corrugated-wall  waveguide  with  a  smooth-walled  wave¬ 
guide  having  the  same  average  radius.  In  the  presence  of 
a  background  plasma,  the  space-charge-limiting  current 
in  a  smooth-walled  cylindrical  waveguide  can  be  in¬ 
creased  by  a  factor  of  ( 1  where  /  =  n  jnt  rep¬ 

resents  the  amount  of  charge  neutralization  provided  by 
the  ions  in  tne  background  plasma.  The  increased  space- 
charge-iimiting  current  may  allow  microwave-generating 
devices  to  transport  more  intense  electron  beams,  leading 
to  the  possibility  of  generating  enhanced  power  levels. 

The  space-charge-  niting  current  in  the  presence  of  the 
plasma  for  a  thin  hollow  beam  with  a  mean  radius  Rb  and 
a  thickness  A  «  Rh  is  given  by  [18]: 


magnetic  field  and  propagating  axially  through  a  slow- 
wave  structure.  The  slow-wave  structure  consists  of  a  cy¬ 
lindrical  waveguide  with  a  periodically  varying  wall  ra¬ 
dius  R(z )  sinusoidally  rippled  about  the  mean  radius  R0, 
such  that: 

R(z)  =  Ro  +  h  cos  koz  (4) 

as  shown  in  Fig.  5,  where  h  is  the  ripple  amplitude,  Zo  is 
the  ripple  period,  and  fco  =  2*/zq. 

The  slow-wave  structure  supports  a  set  of  electromag¬ 
netic  modes  in  the  waveguide  which  has  phase  velocities 
parallel  to  the  beam  velocity  and  propagate  at  less  than 
the  speed  of  light.  These  slow  electromagnetic  waves  in¬ 
teract  resonantly  with  the  negative  energy  slow  space- 
charge  wave  supported  by  the  beam,  which  leads  to  an 
instability  that  transfers  kinetic  energy  from  the  beam  to 


/,  = 


1 7[?2/3  ~  1) 


3/2 


[2  In  (*„/*»)] [1  -/) 


[kA] 
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where  y  =  (1  -  v2/c2 )~ 1  /2  is  the  relativistic  mass  factor 
for  the  beam  electrons.  By  increasing  the  beam  current 
through  introduction  of  a  background  plasma,  a  new  lim¬ 
itation  is  encountered  which  is  imposed  by  the  et  of 
beam-plasma  instability.  This  instability  will  occui  at  cur¬ 
rent  levels  that  are  larger  than  the  vacuum-limiting  current 
by  the  factor  [2]: 


7(1  ~  7~  ) 
(1  -7-2/J). 


3/2 


(6) 


For  mildly  relativistic  electron  beams  having  7  =  3,  the 
actual  current  carried  by  the  beam  in  the  presence  of 
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plasma  may  be  as  much  as  seven  times  larger  than  the 
vacuum  case.  The  increased  injected  current  level  may 
simply  allow  operation  at  higher  beam  power  and  result 
in  enhanced  microwave  power  output  without  affecting  the 
interaction  efficiency  or  the  physical  interaction  mecha¬ 
nism. 

IV.  Linear  Theory  of  Plasma-Loaded  Corrugated 
Waveguides 


field  can  be  expressed  as  [23]: 

E.  =  Z  A„J0  r)  exp  [i{k,,z  -  wf)]  (7) 

"  =  \aq  / 


•  exp  [i(knz  -  w/)]  (8) 


In  addition  to  affecting  the  space-charge-limiting  cur¬ 
rent,  the  presence  of  plasma  in  the  slow-wave  structure 
may  completely  alter  the  nature  of  the  beam-wave  inter¬ 
action,  leading  to  greatly  enhanced  device  power-gener¬ 
ation  efficiency.  The  studies  reported  here  are  the  first  that 
clearly  belong  to  this  latter  category. 

The  linear  theory  of  thin  annular  beams  in  an  evacuated 
corrugated  wall  structure  has  been  treated  both  analyti¬ 
cally  and  numerically  by  researchers  in  the  U.S.  [19]  and 
the  Soviet  Union  [20].  The  introduction  of  plasma  into 
the  slow-wave  structure,  however,  necessitates  the  devel¬ 
opment  of  an  extension  to  the  vacuum  theory.  Previous 
treatment  of  the  linear  theory  of  plasma-loaded  structures 
[21],  [22]  used  the  approximation  of  small  corrugation 
amplitudes  ( h/R0  «  1).  Recently,  a  linear  theory  al¬ 
lowing  for  the  treatment  of  arbitrarily  large  corrugation 
amplitudes  has  been  developed.  A  complete  description 
of  the  theory  can  be  found  in  [23],  The  main  results  of 
this  linear  theory  describing  the  excitation  of  electromag¬ 
netic  waves  in  a  plasma-filled  BWO  by  a  relativistic  elec¬ 
tron  beam  are  presented  below,  where  the  dispersion  re¬ 
lation  and  growth  rates  of  plasma-filled  periodic  structures 
have  been  calculated. 

In  developing  the  linear  theory  the  following  assump¬ 
tions  were  made:  (i)  The  magnitude  of  the  axial-guiding 
magnetic  field  is  taken  to  be  infinite,  confining  the  elec¬ 
tron  motion  in  the  beam  to  one  dimension;  (ii)  the  beam 
is  monoenergetic,  with  an  axial  streaming  velocity  vh;  (iii) 
the  beam  is  of  uniform  density  Nh  and  completely  fills  the 
waveguide;  (iv)  the  waveguide  is  infinite  in  length  and  the 
waveguide  wall  is  perfectly  conducting;  (v)  only  the  sym¬ 
metric  (TM0J  electromagnetic  modes  are  considered;  and 
(vi)  the  plasma-guide  modes  are  neglected. 

In  general,  the  linearized  normal  electromagnetic  modes 
of  a  periodic  corrugated  wall  structure  are  the  transverse 
electric  (TE)  modes,  with  electric  and  magnetic  field 
components  £«,  Br,  and  Bz,  and  the  transverse  magnetic 
(TM)  modes,  with  field  components  £,,  £.,  and  Be.  In  our 
idealized  one-dimensional  system  the  TE  waves  do  not 
perturb  the  axial  electron  motion  and  will  not  be  consid¬ 
ered  further,  since  they  cannot  cause  beam  instabilities. 
The  TM  waves,  on  the  other  hand,  are  capable  of  per¬ 
turbing  the  axial  beam  velocity  and  beam  density  and  as¬ 
sumption  ( v )  above  is  justified. 

The  periodicity  of  the  slow-wave  structure  permits  the 
field  components  and  beam  perturbation  to  be  expanded 
in  an  infinite  Floquet  sum.  Assuming  azimuthal  symme¬ 
try,  the  axial  and  radial  components  of  the  TM0,  electric 


where  wh  is  the  beam-plasma  frequency,  and 


k„  =  k.  +  nkc 


(w  -  k„v)‘ 


Assumption  (iv)  requires  that  the  tangential  electric  field 
vanishes  at  the  waveguide  wall.  Using  the  axial  and  radial 
electric-field  components  of  (7)  and  (8).  we  express  this 
boundary  condition  in  matrix  form  as: 

oo 

D  A  =  Z  A„Dm  „  =  0  (9) 

m.n  =  —  oo 

where  A  is  a  column  vector  with  elements  A„,  and  D  is  a 
matrix  with  elements  Dm  A  dispersion  relation  results 
when  a  nontrivial  solution  to  (9)  is  found  by  solving  the 
determinant  equation  det  |D|  =0. 

In  order  to  numerically  analyze  experimental  systems, 
the  infinite  matrix  of  (9)  must  be  truncated  to  a  finite  size. 
For  example,  to  analyze  an  experimental  system  with  a 
corrugation  amplitude  ratio  h/R0  =  0.3,  a  truncated  ma¬ 
trix  with  a  rank  of  nine  was  found  to  be  sufficient. 

Fig.  6  displays  the  calculated  TM0i  mode  dispersion 
curves  in  a  plasma-filled  corrugated  waveguide  for  a  va¬ 
riety  of  plasma  densities.  The  beam  space-charge  line,  as¬ 
suming  infinitesimal  beam  density  (uh  =  0).  is  superim¬ 
posed  on  this  dispersion  curve  and  has  been  plotted  for 
the  case  which  corresponds  to  the  experimental  parame¬ 
ters  to  be  described  in  the  next  section.  As  can  be  seen  in 
Fig.  6,  the  presence  of  the  plasma  tends  to  raise  the  TMot 
mode  cutoff  frequency  relative  to  the  vacuum  cutoff  fre¬ 
quency  while  also  causing  a  decrease  in  the  group  velocity 
everywhere,  resulting  in  a  “flattening”  of  the  dispersion 
relation. 

When  the  beam  density  is  neglected  (w*  =  0),  the  ra¬ 
diation  frequency  w  is  real  for  real  values  of  the  axial  wave 
number  k.  and  no  instabilities  are  expected.  To  obtain  the 
TM0,  growth  rates,  a  nonzero  value  of  must  be  used 
and  the  dispersion  equation  must  be  solved  for  complex 
values  of  w  for  real  values  of  k..  Fig.  7  displays  the  growth 
rate  Im(w/2ir)  calculated  for  the  experimental  parame¬ 
ters.  The  results  are  similar  to  those  found  in  [15],  with 
the  exception  that  in  this  case  the  peak  growth  rate  is  a 
function  of  the  plasma  density  as  shown  in  Fig.  8.  The 
calculations  predict,  therefore,  a  considerable  decrease  in 
the  linear  growth  rate  for  plasma  densities  of  Np  s  101' 


so: 
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WAVENUMBER  kt  (cm'1) 

Fig.  6.  The  calculated  effect  of  varying  the  plasma  density  on  the  corru¬ 
gated  waveguide  dispersion  (/t0  =  1 .445  cm.  h  -  0.445  cm.  Zo  =  1.67 
cm). 


WAVENUMBER  I,  (cm'1) 

Fig.  7.  Calculated  linear  growth  rate  for  the  TM,,  mode  versus  wavenum¬ 
ber  k.  for  the  case  where  Nr  «  2  •  10"  ctrT\  Nh  «  6.3  •  10'°  cm'1 
(solid  beam). 


PLASMA  OENSITY  N,(,lo"em'*) 


Fig.  8.  Calculated  peak  growth  rate  /„(«)/ 2r  versus  the  plasma  density 
fortheTMo,  mode  (solid  beam  of  density  Af*  =  6.3  •  I0'°cm'3). 

V.  Measurement  of  Efficiency  in  a  Plasma- 
Loaded  BWO 

Early  experiments  testing  the  effect  of  plasma  on  the 
operation  of  vacuum  BWO’s  utilized  background  plasmas 
produced  by  electron-beam  impact  ionization  of  a  low- 
pressure  neutral  background  gas  to  demonstrate  enhanced 


RELATIVISTIC  %  MAGNETIC 


Fig.  9.  Schematic  diagram  of  the  plasma-loaded  BWO  experiment.  The 
electron  beam  is  injected  from  the  left  and  the  radiation  is  extracted  at 
right,  where  the  plasma  is  injected. 


microwave  output  power  [7],  [8].  In  recent  experiments 
[9]  an  independently  controllable  Argon  plasma  source 
was  used  to  inject  plasma  directly  into  the  BWO  slow- 
wave  structure,  as  shown  in  Fig.  9. 

A  hollow  relativistic  electron  beam  wt..<  an  injection 
energy  of  630  keV,  a  beam  current  of  2.3  kA,  and  a  pulse 
duration  of  100  ns  (Fig.  10)  was  propagated  in  a  sinu¬ 
soidally  rippled  slow-wave  structure  immersed  in  a  uni¬ 
form  axial  magnetic  field  of  - 12  kG.  The  slow- wave 
structure  had  an  average  radius  R0  —  1.445  cm,  a  corru¬ 
gation  amplitude  h  ~  0.44  cm,  a  corrugation  period  of  Zo 
=  1 .67  cm,  and  a  structure  length  of  8  periods.  The  elec¬ 
tron  beam  had  an  average  radius  of  0.8  cm,  a  thickness 
of  0.2  cm,  and  an  electron  density  of  -5  x  10M  cm~3. 

A  coaxial  plasma  gun  [24]  was  located  approximately 
100-cm  downstream  of  the  slow-wave  structure  in  a  field 
free  region  and  generated  an  Argon  plasmoid  which 
crossed  the  magnetic  field  lines  at  an  average  velocity  of 
about  1.2  cm /ns  on  its  way  towards  the  interaction  re¬ 
gion.  The  system  was  pre-evacuated  to  a  pressure  of  <4 
x  10~5  T. 

In  separate  experime'  the  plasma  column  density  was 
measured  with  both  a  jj-  and  70-GHz  microwave  inter¬ 
ferometer  and  its  velocity,  density,  and  temperature  was 
measured  with  Langmuir  probes.  Preliminary  results  in¬ 
dicate  that  the  average  plasma  density  could  be  varied  by 
changing  t1’  ■  plasma  gun  voltage  and  gas  pressure.  Real¬ 
izable  plas.  -  densities  in  the  slow-wave  str  ture  ranged 
from  0  up  to  a  maximum  value  estimated  at  -  1012  cm  3. 

There  were  three  principal  overlapping  time  scales  in 
the  BWO  experiments:  The  longest  time  scale  was  on  the 
order  of  10  ms,  corresponding  to  the  quarter-cycle  time  of 
the  pulsed  axial-guiding  magnetic  field.  The  second  long¬ 
est  time  scale  was  on  the  order  of  100  ps,  rresponding 
to  the  plasma  generation  and  propagation  time  from  the 
plasma  gun  to  the  slow-wave  structure.  Finally,  as  noted 
above,  the  beam  e  duration  was  100  ns.  The  three 
time  scales  were  all  synchronized  such  that  the  beam  and 
plasma  interacted  in  the  slow-wave  structure  during  the 
time  that  the  applied  magnetic  field  was  nearly  constant 
and  at  its  peak  value. 

The  plasma  uensity  in  the  slow-wave  structure  was  in¬ 
tentionally  varied  from  shot-to-shot  in  order  to  determine 
its  effect  on  the  operational  efficiency  of  the  BWO.  The 
plasma  density  was  varied  in  the  slow-wave  structure  by 
delaying  the  time  of  beam  injection  into  the  slow-wave 
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Fig.  10.  Diode  voltage  (top  trace)  and  BWO  beam-current  waveforms  used 
in  the  experiment.  An  annular  beam  with  an  average  radius  of  0.8  cm 
and  a  thickness  of  0.2  cm  was  used. 


BEAM  INJECTION  TIME  DELAY  t^sl 

Fig.  1 1.  Peak  microwave  power  (at  8.4  GHz)  versus  the  beam-injection 
time  delay.  By  adjusting  this  delay  we  can  set  the  desired  plasma  density 
in  the  BWO  before  the  injection  of  the  electron  beam. 


structure  relative  to  the  time  of  initial  plasma  generation. 
By  changing  the  relative  delay  times  from  shot-to-shot,  a 
greater  or  lesser  density  plasma  was  allowed  to  diffuse 
into  the  slow  wave  structure  at  the  instant  of  beam  injec¬ 
tion. 

An  enhancement  of  the  microwave  power-generation  ef¬ 
ficiency  was  observed  over  a  wide  range  of  injected 
plasma  densities,  0  <  Np  <  Np(cr),  where  NpUr)  is  a  crit¬ 
ical  plasma  density.  This  critical  plasma  density,  to  be 
discussed  in  Section  VII,  occurred  at  a  delay  time  Ar  - 
90  fts  relative  to  the  initial  generation  of  the  plasma.  The 
efficiency  enhancement  factor  over  vacuum  BWO  opera¬ 
tion  was  found  to  be  dependent  on  the  plasma  density  and 
reached  a  maximum  value  of  eight  when  the  electron  beam 
was  injected  into  an  optimized  plasma  density. 

The  points  of  maximum  BWO  efficiency  enhancement 
were  found  to  occur  for  two  Ji.Tcient  time  delays— once 
for  a  short  time  delay  At  -  60  /as,  which  occurred  during 
the  plasma  buildup  in  the  slow- wave  structure,  and  then 
again  for  a  longer  delay  A/  ~  100  /xs,  which  occurred 
during  the  plasma  decay.  Fig.  11  is  a  plot  of  the  peak 
radiated  power  (at  8.4  GHz)  as  a  function  of  the  beam- 
injection  time  delay  At  gathered  over  a  number  of  shots. 
At  the  points  of  maximum  enhancement,  indicated  by  the 
two  relative  maxima  of  Fig.  1 1,  the  interaction  efficiency 
increased  to  almost  40%,  compared  with  approximately 
5%  for  a  vacuum  BWO  under  the  same  operating  condi¬ 
tions.  The  output  power  was  measured  using  a  side  cou¬ 
pler  calibrated  for  the  TM0,  mode  over  the  relevant  band¬ 
width. 

The  plasma  density  in  the  slow-wave  structure  reached 
its  highest  values  for  beam-injection  time  delays  in  the 
range  60  ns  <  At  <  100  n s.  In  this  range  the  BWO  ra¬ 
diation  at  8.4  GHz  from  the  fundamental  TM0)  mode  was 
quenched  and  higher  frequency  emission  in  the  12-18- 
GHz  range  was  observed.  At  these  higher  plasma  densi¬ 
ties  the  BWO  was  considered  to  be  overdriven  and 
switched  from  the  fundamental  TM0|  mode  to  TM02  and 
possibly  higher  order  modes,  as  indicated  in  Fig.  12. 

An  example  of  the  change  in  frequency  spectra  that  ac¬ 
companies  mode  switching  is  shown  in  Fig.  13,  which 
plots  the  single-shot  experimental  spectrum  for  an  effi¬ 
ciency-enhanced  BWO,  shown  in  Fig.  13(a),  and  an  over¬ 
driven  BWO,  as  shown  in  Fig.  13(b).  These  experimental 
curves  were  obtained  using  a  dispersive  line  technique, 
where  different  frequency  components  are  resolved  by 
their  different  propagation  velocities  in  a  dispersive  wave- 
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Fig.  12.  Microwave  radiation  detected  from  the  plasma-loaded  BWO  in 
the  frequency  bands:  (a)  12-18  GHz  (corresponding  to  emission  at  high- 
order  mode);  and  (b)  18-26  GHz  (which  may  correspond  to  emission 
due  to  an  eiectromagnetically  pumped  FEL;  see  discussion). 
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Fig.  13.  Spectral  results  from  plasma-enhanced  backward- wave  oscillator, 
(a)  Maximum  enhancement;  and  (b)  overdriven  device.  The  "prompt 
signal"  is  a  timing  reference. 


guide.  Fig.  13(a)  plots  the  power  spectrum  emitted  by  the 
BWO  under  the  condition  of  strong  plasma  enhancement, 
where  the  dominant  frequency  component  is  centered  near 
8.4  GHz.  A  timing-reference  microwave  signal  (labeled 
prompt)  is  superimposed  on  the  dispersed  signal.  Fig. 
13(b)  represents  the  spectrum  for  an  overdriven  device  in 
which  mode  switching  has  clearly  occurred  and  frequency 
components  in  the  range  of  12-18  Gz  are  present,  corre¬ 
sponding  to  TM02  or  higher  order  modes  as  well  as  plasma 
oscillation.  At  even  higher  plasma  densities  microwave 
breakdown  within  the  device  will  occur. 

Strong  microwave  emission  was  also  detected  in  the 
18-26  GHz  band,  proportional  in  amplitude  to  the  fun- 
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damental  TM0I  backward-wave  oscillations.  These  high 
frequencies  are  believed  to  be  produced  by  a  free-electron 
laser  interaction  driven  by  an  electromagnetic  pump,  as 
previously  reported  125},  [26]. 

In  our  experiments  the  introduc'ion  of  plasma  into  the 
device  had  no  adverse  effects  on  the  diode  voltage  and 
current.  The  beam  current  entering  the  slow- wave  struc¬ 
ture  was  unaffected  by  the  presence  of  the  plasma  and  no 
effect  on  the  diode  shorting  time  was  observed. 

VI.  Discussion  of  Results,  Coupling  Mechanisms, 
and  Future  Trends  in  Plasma  Microwave 
Electronics 

As  the  experiments  described  in  the  previous  section 
demonstrate,  the  presence  of  a  controlled  plasma  in  the 
BWO  slow-wave  structure  strongly  enhances  the  interac¬ 
tion  efficiency  over  a  very  wide  range  of  plasma  densities . 
In  addition  to  efficiency  enhancement,  two  other  interest¬ 
ing  features  of  the  experimental  results  are:  BWO  plasma- 
saturation  effects,  and  mode  switching. 

It  vas  anticipated  that  overdriving  the  device  with 
plasma  above  a  critical  density  Np(er)  would  lead  to  a 
quenching  of  the  BWO  interaction.  As  discussed  in  Sec¬ 
tion  IV,  the  presence  of  plasma  in  the  slow-wave  structure 
tends  to  raise  the  lower  cutoff  frequency  while  simulta¬ 
neously  causing  a  flattening  of  the  overall  dispersion 
curves.  If  the  device  is  overdriven  to  a  point  where  the 
lower  and  upper  cutoff  frequencies  are  similar,  the  elec¬ 
tromagnetic-mode  group  velocity  is  drastically  reduced 
and  no  backward  wave  interaction  is  possible. 

The  linear  theory  of  electromagnetic  waves  in  plasma- 
loaded  BWO’s,  presented  in  Section  III,  predicts  that 
BWO  TM0|  quenching  by  plasma  saturation  might  be  ob¬ 
served  for  critical  densities, 

%(fr)  s  10l2[cm'3].  (10) 

If  the  plasma  density  is  above  this  critical  value,  the  linear 
gr  wth  rate  will  begin  to  decrease  as  shown  in  Fig.  8. 
Linear  theory  also  predicts  that  mode  switching  may  oc¬ 
cur  at  about  the  same  critical  density  [23],  Indeed,  both 
BWO  plasma  saturation,  as  seen  in  Fig.  11,  and  mode 
switching,  as  seen  in  Figs.  12  and  13,  were  observed  when 
the  device  was  overdriven.  This  occurred  when  the  beam- 
injection  delay  time  was  adjusted  to  At  -  90  pis,  corre¬ 
sponding  to  an  estimated  plasma  density  of  Np  <  2  x  101 1 
cm'3,  which  is  far  below  the  value  predicted  by  the  linear 
theory. 

The  linear  theory  for  a  plasma-loaded  fundamental 
TM0i  mode  in  a  BW  predicts  that  the  interaction  fre¬ 
quency  of  the  BWO  ill  increase  with  increased  p  isma 
density,  as  may  be  seen  in  Fig.  6.  This  predicted  increase 
has  not  yet  been  observed  experimentally.  Recent  data  in¬ 
dicates  that  the  radiation  frequency  remains  fixed  (within 
the  measurement  resolution)  despite  changes  in  the  plasma 
density.  The  growth  rate  predicted  by  the  theory  i?  only 
slightly  larger  in  the  presence  of  the  plasma  as  cor  red 
to  the  vacuum  BWO  case. 

Thus  while  the  linear  theory  is  successful  at  predicting 
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some  experimentally  observed  trends  such  as  reduction  of 
the  growth  rate  (plasma  saturation)  and  mode  switching 
due  to  the  presence  of  the  plasma,  it  does  not  produce 
numerical  results  that  are  reasonably  close  to  the  experi¬ 
mental  data.  This  leads  us  to  the  conclusion  that  the  pres¬ 
ent  linear  theory  must  be  extended  to  account  for  the  ob¬ 
served  experimental  results. 

One  physical  mechanism  which  may  explain  the  en¬ 
hanced  efficiency  is  a  three-wave  interaction  involvir^  the 
beam,  electromagnetic,  and  T-G  plasma  waves  in  the  pe¬ 
riodic  slow-wave  structure.  In  a  nonperiodic  structure,  the 
normal  T-G  modes  will  not  interact  with  a  relativistic 
electron  beam,  as  their  phase  velocities  are  small  com¬ 
pared  to  the  phase  velocity  of  the  beam  space-ch  rge 
wave.  In  a  periodic  structure,  however,  the  T-G  modes 
are  also  expected  to  exhibit  periodicity  in  k-space,  as  re¬ 
quired  by  the  Floquet  theorem.  The  dispersion  curve  for 
a  periodic  plasma  is  quite  different  from  that  shown  in 
Fig.  1(b)  for  a  plasma  in  a  smooth-walled  waveguide. 
Most  significantly,  there  exist  some  periodic  T-G  modes 
that  have  phase  velocities  that  are  equal  to  that  of  the  re¬ 
lativistic  electron  beam  and  are  capable  of  supporting 
backward  waves. 

As  an  example.  Fig.  14  plots  the  uncoupled,  shifted 
dispersion  curves  for  the  T-G(0  2)  modes  (labeled  as 
“piasma  waves”),  and  also  the  approximate  coupled  dis¬ 
persion  relation  for  this  mode.  The  dispersion  curve  for 
the  lowest-order electromagnetic  mode  (labeled  ”TM0i”) 
can  also  be  seen.  Azimuthal  symmetry  has  been  assumed, 
and,  for  clarity,  only  the  first  five  periods  in  k  have  been 
plotted. 

For  the  range  of  values  of  the  plasma  density  and  beam 
energy  relevant  to  the  experiments  reported  in  this  paper, 
the  beam  is  simultaneously  in  synchronism  with  the  back¬ 
ward  branch  of  the  electromagnetic  structure  wave  and 
several  backward  branches  of  the  periodic  plasma  waves. 
The  slow  space-charge  wave  and  the  lowest  order  TM0! 
electromagnetic  wave  have  the  same  phase  velo.  ,ty  u/k 
as  do  one  or  more  of  the  periodic  ~-G  waves.  In  addition, 
the  three  waves  all  share  a  stron..  common  componer  ff 
the  axial  electric  field  £..  It  is  possible  that  the  enhanced 
efficiency  observed  in  the  experiment  results  from  this 
three-wave  synchronism;  namely,  induced  scattering  of 
the  electromagnetic  radiation  of  electrons  in  an  electro¬ 
static  field  produced  by  the  background  plasma  in  the 
presence  of  a  periodic,  perfectly  conducting  wall. 

An  efficient  three-wave  interaction  also  implies  a  very 
broad  gain  curve  with  respect  to  the  electron-beam  en¬ 
ergy.  This  implication  is  consistent  with  the  foregoing 
proposed  mechanism,  since  the  beam  electrons  tend  to 
stay  in  synchronism  with  both  of  the  backward  branches, 
even  as  the  beam  is  losing  energy.  A  change  of  40%  in 
the  beam  kinetic  energy  corresponds  to  a  change  of  less 
than  5%  in  the  streaming  velocity  for  our  experimental 
conditions. 

A  second  mechanism  which  may  explain  the  enhanced 
emission  in  the  presence  of  the  plasma  is  a  stimulated  Ra¬ 
man  scattering  process  which  may  be  described  using  a 
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Fig.  14.  Schematic  sketch  of  the  dispersion  diagram  for  the  electromag¬ 
netic  TMoi  mode,  the  plasma-guide  T-G(0  , ,  mode,  and  a  space-charge 
wave  on  a  low-density  beam  traveling  at  vh  =  0.9c  (see  discussion). 


Fig.  15.  Schematic  sketch  of  Stokes  diagram  for  three-wave  interaction, 
which  involves  a  backward  branch  of  the  plasma-guide  (T-G)  mode,  and 
two  TMo,  waves  of  somewhat  different  frequencies. 


Stokes  diagram,  as  shown  in  Fig.  15.  In  this  figure,  the 
electromagnetic  wave  (point  A)  decays  into  a  backward 
plasma-guide  (T-G)  wave  and  an  electromagnetic  wave 
which  is  slightly  shifted  in  frequency.  Thus  this  mecha¬ 
nism  proposes  a  multiwave  interaction  as  the  explanation 
for  enhanced  efficiency. 

It  should  be  noted  that  for  the  experiments  described  in 
this  work,  the  electron-beam  current  was  always  below 
that  of  the  vacuum  limit  for  the  device  geometry.  As  a 
result,  we  have  not  yet  studied  the  radiation  generation 
mechanism  at  beam  currents  above  this  limit.  As  dis¬ 
cussed  in  Section  III.  the  introduction  of  plasma  into  the 
slow-wave  structure  may  have  the  effect  of  neutralizing 
the  RF  space  charge  in  the  beam  (allowing  a  greater  de¬ 
gree  of  axial  bunching)  and  also  may  allow  the  propaga¬ 
tion  of  beam  currents  above  the  vacuum  space-charge 
limit.  This  phenomenon  is  of  great  interest  to  the  field  of 
relativistic  microwave  electronics  and  may  be  the  subject 
of  future  experimental  work. 

It  is  anticipated  that  plasma  injection  may  also  prove 
beneficial  to  a  variety  of  high-power  microwave  devices 
1 12]— [  16],  [21].  Gyrotrons,  for  example,  may  achieve  a 
greater  degree  of  tunability  and  yield  higher  output  power 
by  overcoming  the  space-charge  limit  [14]. 

It  was  predicted  theoretically  [12]  that  larger  gain  can 


be  realized  in  plasma-loaded  free-electron  lasers  that  are 
expected  to  exhibit  greater  tunability  and  performance  en¬ 
hancement  by  changing  the  refractive  index  n  in  the  cavity 

(13] - 
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Abstract— A  linear  theory  of  the  excitation  of  electromagnetic  waves 
in  a  plasma-filled  corrugated-wall  waveguide  with  an  arbitrarily  large 
sinusoidal  corrugation  has  been  derived  and  analyzed  numerically.  The 
theory  predicts  that,  when  driven  by  an  electron  beam,  the  presence 
of  a  plasma  in  the  slow  wave  structure  will  cause  an  increase  in  the 
oscillation  frequency,  and  that  the  temporal  growth  rates  of  a  high- 
frequency  mode  approach  those  of  fundamental  mode  for  high  plasma 
densities.  The  latter  result  may  account  for  the  high-frequency  modes 
observed  in  our  plasma-filled  backward-wave  oscillator. 


I.  Introduction 

THE  ELECTRON  beam  excitation  of  slow  wave  struc¬ 
tures  has  been  an  active  research  area  for  more  than 
halt  a  century  (I).  With  recent  progress  in  intense  relativ¬ 
istic  electron  beam  (IREB)  technology,  a  new  generation 
of  high-power  microwave  sources  is  being  developed.  The 
backward-wave  oscillator  (BWO)  is  one  of  the  most 
prominent  IREB-driven  high-power  microwave  sources 
[2)-[51.  Peak  powers  in  excess  of  3  GW  in  a  single-stage 
device,  15  GW  in  a  two-stage  device,  long  pulse  lengths 
of  several  microseconds,  high  efficiencies  (  >50%  ).  and 
a  pulse  repetition  rate  of  400  Hz  have  been  reached  in  a 
variety  of  experiments  (6).  More  commonly,  power  levels 
of  hundreds  of  megawatts  at  frequencies  ranging  from  8 
to  35  GHz  have  been  reported,  although  frequencies  as 
high  as  150  GHz  have  been  demonstrated 
Recent  experiments  at  the  University  of  Maryland  have 
demonstrated  high-efficiency  high-power  operation 
through  the  introduction  of  a  background  plasma  into  an 
X-band  BWO  structure  [7),  (8). 
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As  a  first  step  toward  understanding  this  phenomenon, 
a  linear  dispersion  relation  for  a  plasma-loaded  BWO  has 
been  derived  and  analyzed  numerically,  with  special  em¬ 
phasis  placed  on  using  system  parameters  consistent  with 
experiments  performed  at  the  University  of  Maryland  (7). 
[8].  The  development  of  the  dispersion  relation  may  lead 
to  a  discussion  of  a  three-wave  interaction  (beam/plasma/ 
slow  electromagnetic  wave)  as  a  possible  physical  mech¬ 
anism  to  explain  experimental  observations. 

Section  II  of  this  paper  contains  the  ’  ■  i  vation  of  a  dis¬ 
persion  relation  for  waves  in  an  infinitely  long  plasma- 
loaded  conrugated-wall  waveguide  excited  by  an  electron 
beam.  Section  III  presents  the  numerical  analyses  of  the 
interaction,  and  a  discussion  of  this  work,  along  with  con¬ 
clusions.  is  presented  in  Section  IV. 

II.  Dispersion  Relation  for  a  Plasma-Loaded 
Corrugated-Wall  Waveguide 

The  analyses  presented  here  are  based  on  the  slow-wave 
structure  shown  in  Fig.  I.  The  slow-wave  structure  con¬ 
sists  of  an  axially  symmetric,  cylindrical  waveguide 
whose  wall  radius,  R(z).  varies  sinusoidally  according  to 
the  relation 


R(z)  =  Ro  +  h  cos  ( k0z ) 


(1) 


where  h  is  the  corrugation  amplitude.  k{]  =  2ir/c(,  is  the 
corrugation  wavenumber,  and  z0  is  the  length  of  the  cor¬ 
rugation  period. 

A  solid  uniform  electron  beam  of  density  Nh  and  radius 
Rh  <  R0  -  h  goes  through  the  waveguide,  which  is  loaded 
completely  with  a  cold,  uniform,  collisionless  plasma  of 
density  Nr.  Th**  entire  system  is  immersed  in  a  strong, 
longitudinal  magnetic  field,  which  magnetizes  both  the 
beam  and  the  plasma.  There  is  no  restriction  on  the  rela¬ 
tionship  between  the  beam  density  and  the  plasma  den¬ 
sity, 
quency. 

respectively,  are  assumed  to  be  much  snuller  than  the 
electron  cyclotron  frequency. 

We  further  assume  that  i)  the  beam  current  is  taken  to 
be  much  less  than  the  space-charge-limiting  current  for  a 
smooth-walled  waveguide  with  the  same  average  radius: 
ii)  the  beam  is  monoenergetic  with  all  electrons  having 
identical  longitudinal  velocity  r:  and  iii)  the  beam  has  an 


sity.  but  the  plasma  frequency  and  the  beam  plasma  fre¬ 
quency.  =  (e''Np/m,.t o)^'  and  uih  =  (e: Nh/mft0)'". 
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CORRUGATED  WALL 

WAVEGUIDE 

R(r)=  R0*  ncosk0z- 


■  UNIFORM  PLASMA  N„ 


r SOLID  BEAM 
I  N„ 


1  0  0 

[(}  =  0  10 
0  0 


and  by  [10] 


73(oj  -  k„v)' 


where  7  is  the  relativistic  factor,  and  k„  is  the  wavenum¬ 
ber  (to  be  defined  later). 

The  periodicity  of  the  corrugated  structure  permits  the 
field  components  and  the  beam  perturbation  to  include  an 
infinite  number  of  components  in  the  wavenumber  space 
according  to  the  Floquet  theorem  [11],  Only  TM  waves 
perturb  the  axial  velocity  and  the  density  of  the  beam. 
Accordingly,  we  will  focus  our  attention  on  the  dispersion 
relations  for  the  axisymmetric  case  of  these  modes  alone. 
From  Maxwell’s  equations  and  (2)  and  (3),  one  van  obtain 
the  following  expressions  for  the  axial  and  radial  RF  elec¬ 
tric  fields  £.  and  £,  for  the  axisymmetric  TM  electromag¬ 
netic  modes: 


£,(r.  z)  =  L  £.„  exp  [/(*,,;  -  ut) J 


_  v 

-  A— >  1 

n  ~  -  oo  (j)  “  ^ 

C~ 


exp  [i(k„z  —  cur ) 


where  k.  is  the  axial  wavenumber,  and 
kn  =  k.  +  nko 


V2"/ko 

Fig  I  Slow  wave  slruciure  and  beam  model  of  secnons  1I-1V  lor  solid 
electron  beam  of  radius  R,,  within  perfectly  conducting  plasma-loaded 
corrugated-wall  waveguide. 


equilibrium  position,  and  is  free  from  instabilities  which 
cause  macroscopic  deformations. 

The  dispersion  relation  for  a  plasma-filled  BWO  was 
previously  derived  for  the  case  of  a  slow-wave  structure 
with  a  small  comigation  amplitude  [9],  (10)  and  only  in¬ 
cluded  terms  up  to  the  order  of  (h/R0)'-.  In  our  experi¬ 
ments,  however,  the  corrugation  amplitude  was  not  small 
{h/R0  =  0.308),  necessitating  the  derivation  of  a  disper¬ 
sion  relation  for  arbitrarily  large  values  of  h /R0  <  1 . 

The  assumption  of  an  infinitely  large  guiding  magnetic 
field  implies  that  the  only  nonvanishing  terms  in  the  rel¬ 
ative  dielectric  tensor  are  along  the  major  diagonal.  For 
the  beam-plasma  case  in  a  linearized  scheme  of  treat¬ 
ment,  the  relative  dielectric  tensor,  in  cylindrical  coordi¬ 
nates,  may  be  given  by 


=  A, Vo  (^r)  0  g  r  g  Rh 


Rh  §  r  ^  R(z) 


Xn  —  Rq  (  2  ^  w  )  1  ~  ~  — - - 3 

\c  /  w  7  (w  ~  £,.*’) 


yr,  -  *o  I  - 


1  -  -  ■ 

U1  } 


The  constants  B„  and  C.  in  (6)  can  be  expressed  in  terms 
of  the  constant  A„  fn  the  conditions  for  the  continuity 
of  £,  and  £-  ar  r  -  R  .» 

R»  ^  KnAn%  Cn  ^  Ln An 

K„  =  ynV0(xn6)A/,(yn6)  -  x„J,  (xn6)N0(  y„6)  'j 

E„  =  x„J0{y„6)J^xnb)  -  y„J0(x„6)  •/,(>', ,$)j 

6  =  Rh/R o- 

The  RF  electric  fields  £.  and  Er  must  satisfy  the  bound¬ 
ary  condition  that,  at  the  perfectly  conducting  corrugated 
waveguide  surface,  the  tangential  electric  field  £,  must  be 
zero;  i.e., 

E,{r  =  J?(z)) 

«  £-( r  =  R(z))  +  Er(r  =  R(z))  =  °- 

az 


Substituting  (4)-(9)  into  (10),  the  dispersion  relation 
which  describes  the  dependence  of  the  angular  oscillation 
frequency  «  on  the  axial  wavenumber  k.  is  derived  from 
the  requirement  for  the  existence  of  a  nontrivial  set  of 
solutions  of  amplitudes  A„. 

Equation  (10)  is  a  dispersion  relation  which  gives  a  re¬ 
lationship  between  u>  and  k.\  however.  (10)  is  only  of  lim¬ 
ited  use  in  itself,  because  it  is  derendent  on  the  coordinate 
;  and  involves  as  yet  unknown  oefficients  A„.  To  elimi¬ 
nate  these  quantities,  we  proceeo  as  in  [9]  and  [10],  mul¬ 
tiplying  (10)  by  exp  (  -imk^z)  and  integrating  from  z  = 
-ir /k0  to  ;  =  ir /kn.  The  important  difference  from  pre¬ 
vious  works  [9],  [10]  is  that  this  work  uses  an  arbitrary 
comigation  amplitude  h,  and  all  the  higher  order  terms 
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TABLE  I 

Numerical  Factors  W(  I  h  -  m  I.  v  i  as  Defined  by  Eol  a  rios  <  1 5 >  in  Ti  \  r" 


\s 

0 

LJ 

2 

3 

4  X 9 

6 

A 

7 

k 

10 

0 

1  ^ 

r 

4  v 

r 

2304^ 

r - 

147436^ 

r 

14749600 

'  (j 

2 

/ 

16 

364 

18432 

1474660 

2  > 

r  ^ 

a 

/ 

98 

3072 

184320 

17694720 

3 

Y - 

/j 

Lx"1 

4a 

768 

30720 

2211840 

4 

7^ 

1  r 

k 

364 

7680 

366640 

30965760 

5 

3840 

92160 

3180980 

6 

1 

?  J 

k 

aaoao 

1290240 

82373360 

7 

V. 

M 

k 

649120 

20643640 

a 

_ 

_ 

_ 

CO 

i _ 

10321920 

371589120 

Blank  boxes  indicate  lhai  I  /Hi  i  n  -  in  |.  s )  =  0  Factors  *'  t  3  x  3.  5  x  5.  and  9  x  9  deierminani  calculations  arc  shown 
by  thick  lines  Here  s  is  given  by  s  =  2q  +  |  u  -  m  |  as  in  equation  ( 1 5).  Adjacent  diagonal  numbers  are  related  lo  each  other 
by  a  factor  enclosed  by  a  circle 


are  included  in  the  analysis.  The  resulting  dispersion  re¬ 
lation  is  found  to  be 


x 


+ 


dz  =  0. 


(ID 

To  solve  (1 1),  J0(ynR(z)/R0)  and  N0(  ynR(z) / Ro)  are 
Taylor  expanded  about  R(z)  =  fly.  The  following  rela¬ 
tions  are  useful  formulas  which  facilitate  the  calculation: 


cos 


25 


1  /2.v 


‘“2^[,h{j)COs2h-J,“  +  2  \s 


2\  *  I  I  V1 

cos  U  =  —  2a 

2*  y  =  0 


25  +  I 


J 


COS  [  2  (  5  —  j)  +  I  ]  u 


L 

i:.- 


cos  mu  cos  nu  du  = 


0  (m  *  n) 
ir  (m  =  n  *  0) 


sin  mu  cos  nu  du  =  0. 


Equation  (11)  may  be  expressed  in  matrix  notation: 

ao 

D  A  =  E  D,„„A„  =  0  (12) 

in  ■  n  —  -  ao 

where  A  is  a  column  vector  with  elements  A„,  and  D  is  a 
matrix  with  elements  D„„  obtained  from  the  integral  in 
(II). 

Dmn  =  [  1  +  (n  -  m)Q„](K„C Jm„  +  Z.„Ci)  ( 13) 


where 


rj  _  y  \y»a 

^  inn  —  ^  ” 


P  ' 


'(>■„) 


9  =  0  2 2,<*u,-""q\{q  +  | n  -  mj)! 


Qn  ~  2 


u 


kl 


h 

a t  —  — 

Ro 


(14) 


and  C*„  is  obtained  by  replacing  J0  in  CJ„„,  by  N0. 

Equation  (12)  represents  an  infinite  set  of  linear  equa¬ 
tions  in  unknown  /4,,’s.  Each  term  CJnms  in  (14)  is  of  the 
form: 

r,  _  (y„a)Vo"(.v„) 

H(\n  -  m  | ,  5 ) 


where  s  =  2q  +  |  n  —  m  |  =  1.2.  •  •  •  is  an  integer. 
/0V,(  y„ )  is  the  5  th  derivative  of  the  Bessel  function  y0(  y„ ). 
The  numerical  factor  H ( |  n  -  m  | .  5 )  has  been  calculated 
and  tabulated  in  Table  I.  For  the  diagonal  terms.  H(\n  - 
m|.  5)  in  the  matrix  Dare  found  by  using  (n  -  m)  =  0 
in  the  table.  Similarly,  for  the  first  off-diagonal  terms,  (n 
-  m  =  ±  1  )  is  used,  and  so  on. 

The  dispersion  relation  for  the  plasma-filled  BWO  is 
expressed  by  the  determinant  equation  D  =  det  [D]  =  0. 
In  the  limit  of  infinitesimally  small  /t.  the  dispersion  re¬ 
lation  given  by  (13)  agrees  with  previous  calculations  (9], 
{ 10).  Beam  and  plasma  effects  are  included  in  ( 13)  through 
the  terms  x„  and  y„  given  by  (7)  and  (8),  which  are  the 
functions  of  and  u>h.  For  the  case  without  plasma  (ur 
=  0),  the  calculated  dispersion  curves  coincide  with  those 
previously  calculated  for  the  vacuum  BWO's  [  1 2)— [  15). 

When  no  beam  is  present  (uh  =  0).  Im(w)  becomes 
zero  for  real  values  of  k.,  as  there  can  be  no  instability, 
and  the  numerical  computation  of  (13)  is  simplified.  In 
this  case.  Re(uj).  at  the  intersection  of  the  dispersion  re¬ 
lation  of  the  slow-wave  structure  and  the  beam  space- 
charge  wave  (BSCW),  is  considered  to  be  the  eigenfre- 
quency  of  the  plasma-filled  BWO. 
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As  is  well  known,  waves  with  frequencies  less  than  the 
plasma  frequency  air  can  be  propagated  in  a  smooth- 
walled  cylindrical  plasma  waveguide  ( 1 6j .  This  wave  is 
called  the  Trivelpiece-Gould  mode.  In  the  following  nu¬ 
merical  analyses  of  our  dispersion  relation,  we  coniine 
ourselves  to  electromagnetic  modes  with  angular  free  n- 
cies  greater  than  the  plasma  frequency  un. 

III.  Numerical  Results 

A.  Dispersion  Relation  for  Infinitesimal  Beam  Densities 

First,  we  present  numerical  computations  for  the  case 
of  infinitesimally  small  beam  densities  (o>„  =  0)  where 
the  oscillation  frequency  a>/2x  is  real  for  a  real-valued 
wavenumber  k:.  In  this  case,  y„  =  x„  from  (7)  and  (8). 
and  hence.  K„  =  -2/ttb  and  L„  =  0  in  (9).  The  terms 
including  L„  are  necessary  to  take  into  account  when  the 
beam  density  is  high.  The  following  parameters,  corre¬ 
sponding  to  experimental  conditions,  were  used  in  the 
computations:  relativistic  factory  =  2.23  (i.e..  v/c  = 
0.9):  average  waveguide  radius  R0  =  1 .445  cm;  and  axial 
comr.-  ition  period  20  =  1.67  cm.  The  corrugation  ampli¬ 
tude  was  varied  over  a  wide  range,  and  included  the  ex¬ 
perimental  value,  h  =  0.445  cm.  The  plasma  density  Nr 
was  treated  as  a  variable. 

In  Fig.  2.  the  dispersion  curves  calculated  using  a  3  x 
3 (  —  I  <  [m,  n\  <  1 ),  5  x  5 ( —  2  <  \m.  n)  <  2).  and 
9  x  9(-4  <  [m.  n]  <  4)  determinant,  as  defined  by 
(13),  are  plotted  over  a  range  of  wave  numbers,  k.  =  0 
to  k.  =  7.5248  cm'1,  which  corresponds  to  two  full 
wavenumber  periods  of  the  slow-wave  structure.  The 
plasma  density  was  Nr  =  2  x  10"  cm'1  and  the  corru¬ 
gation  amplitude  was  h  =  0.445  cm. 

For  each  of  the  three  cases  (3  x  3.5  x  5,  and  9x9 
matrices)  shown,  respectively,  by  dashed,  dot-dashed, 
and  solid  lines,  the  terms  with  constant  factors  used  in 
CJmm  in  (15)  of  the  determinant  equation  are  highlighted 
by  the  thick  lines  in  Table  I.  The  three  lowest  curves  in 
Fig.  2  are  labeled  the  ”TM0,  mode”  because  they  are 
related  to  their  well-known  lowest-order  counterpart  in  the 
st  th-walled  plasma-filled  waveguide.  On  the  other 
har.a,  the  next  higher  frequency  mode  to  the  TM0i  mode 
is  labeled  the  ”A”  mod"  which  seems  to  be  composed 
by  displaced  TM01  modes  in  a  smooth-wailed  waveguide. 
The  third  high-frequency  mode  is  named  the  ”B”  mode 
in  the  present  paper.  Other  high-frequency  modes  were 
not  calculated  for  presentation  in  this  paper. 

The  beam  space-charge  wave  (BSCW)  v/c  =  0.9  is 
superimposed  on  the  electromagnetic  modes  in  Fig.  2  and 
is  indicated  by  a  straight  line.  The  intersections  of  the 
BSCW  with  the  TM0i,  A  and  B  electromagnetic  modes 
are  considered  to  be  the  eigenfrequencies  of  the  plasma- 
loaded  corrugated  waveguide. 

In  the  case  of  the  3  x  3  determinant  equation,  only  the 
terms  -  I  <  \m.  n]  <  1  were  used  in  (13).  In  other 
words,  three  displaced  periodic  modes  were  superposed 
and  the  terms  up  to  the  order  of  a;  =  (h/R0):  were  in¬ 
cluded.  It  should  be  noted  that  in  the  case  of  the  3  x  3 
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Fig.  2.  Dispersion  curves  calculated  using  3  x  3.  5  x  5.  and  9x9  de¬ 
terminant  equations  shown,  respectively,  by  dashed,  dot-dashed,  and 
solid  lines  lor  three  electromagnetic  modes.  N,,  =  2  x  10"  cm  ' .  h  = 
0.445  cm.  Straight  line  is  the  beam  space  charge  wave  for  /■/<’  =  0  9. 


determinant,  the  required  periodicity  (from  the  applica¬ 
tion  of  the  Floquet  theorem)  cannot  be  satisfied  even  for 
a  single  wavenumber  period.  0  <  k.  <  3.7624  cm*1,  as 
shown  in  Fig.  2.  This  indicates  that  the  previously  de¬ 
rived  dispersion  relation  |9).  [10]  may  not  be  used  to  ac¬ 
curately  describe  BWO’s  with  large  corrugation  ampli¬ 
tudes  (a  -  I ).  Even  using  a  5  x  5  determinant  equation, 
the  required  periodicity  in  the  dispersion  curve  is  not  ob¬ 
tained  for  the  range  of  wavenumbers  0  <  k-  <  7.5248 
cm'1. 

To  get  satisfactory  periodicity  in  this  wavenumber 
range,  it  has  been  found  that  a  9  x  9  determinant  equa¬ 
tion,  including  terms  up  to  (  y„a)l0./0l0,(  >’«)  in  (14).  must 
be  used.  The  rest  of  the  numerical  results  described  in  this 
paper  were  calculated  using  such  a  9  x  9  determinant. 

The  dispersion  relation  for  extremely  smr  >rrugation 
(h)  values  is  shown  in  Fig.  3  for  the  case  .  Np  =  2  x 
I011  cm*  .  The  dashed  and  the  solid  lines  correspond  to 
h  =  0.01  and  h  =  0. 1  cm.  respectively.  It  can  be  seen 
that  the  TM0i  and  A-modes  are  closely  approximated  by 
the  superposition  of  the  displaced  TM0i  mode  in  a  smooth- 
walled  waveguide,  whereas  the  B-mode  consists  of  dis¬ 
placed  TM0|  and  TM0;  modes  in  a  smooth-  ded  wave¬ 
guide.  The  three  modes  are  increasingly  separated  in  fre¬ 
quency  as  h  increases. 

In  Fig.  4,  the  dispersion  relation  of  the  three  modes  for 
a  various  plasma  densities  Nr,  as  well  as  the  vacuum  case, 
are  shown  for  h  =  0.445  cm.  which  corresponds  to  our 
experimental  parameters  [7],  [8].  The  emr"  waveguide 
results  (dotted  line  in  Fig.  4)  coincide  with  those  in  pre- 
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Fig  3.  Dispersion  curses  lor  N.,  =  2  x  10"  cm  ,V„  =  0  lor  cases  ol  h 
=  0  01  cm  (dashed  line),  and  A  =  0.1  cm  (solid  line). 


h  r  0  445(cm) 
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Fig.  4.  Dispersion  curves  for  three  modes  for  vanous  plasma  densities  iV,, 
In  all  cases.  N,,  *  0  and  h  =  0.445  cm. 

vious  calculations  |12],  (15).  It  is  clear  from  the  figure 
that  the  plasma  presence  in  the  waveguide  shifts  the  dis¬ 
persion  curve  of  the  TM0,  mode  up  in  frequency  also  tends 
to  “flatten''  it;  i.e..  reduces  the  overall  group  velocity  of 
the  curve. 

For  the  A-mode.  the  dispersion  curve  is  very  flat  orig¬ 
inally  ( Nn  =  0 )  and  tends  to  change  with  frequency  as  Nr 
increases.  In  the  same  figure,  the  beam  space-charge  line 
corresponding  to  v/c  =  0.9  is  shown.  The  intersections 
of  this  beam  line  with  the  electromagnetic  modes  are  the 
expected  eigenfrequencies  of  the  plasma  BWO  interaction 
for  infinitesimal  electron  beam  densities. 


The  eigenfrequencies  of  three  modes  versus  plasma 
density  are  plotted  in  Fig.  5.  The  eigenfrequencies  in¬ 
crease  almost  linearly  with  the  plasma  density.  The  dashed 
line  is  the  background  plasma  frequency,  and  it  is  esti¬ 
mated  to  be  approaching  gradually  to  the  eigenfrequency 
of  the  TM0|  mode  for  high  densities,  when  the  resonant 
enhancement  of  the  Raman  scattering  process  was  pre¬ 
dicted  [10]. 

Once  the  dispersion  relation  D  =  0  is  solved  numeri¬ 
cally  (as  in  (13)),  one  can  determine  the  relative  magni¬ 
tudes  between  the  eigenvectors  A„.  where  n  =0.  ±1, 
+  2.  ■  •  •  .  In  other  words,  the  relative  amplitudes. 
\A„/A0\,  of  the  shifted  periodic  modes  can  be  calculated 
for  practical  parameters.  Since  we  have  solved  a  9  x  9 
determinant  equation  to  get  the  dispersion  relation,  in¬ 
homogeneous  simultaneous  equations  of  eight  unknown 
variables  A,J A0  with  n  =  ±1.  •  •  •  .  ±4  can  be  solved 
numerically. 

An  example  of  a  solved  set  of  \A„/A{)\'s  is  shown  in 
Fig.  6  for  the  case  of  Nr  =  2  x  10"  cm~\  From  Flo- 
quet's  theorem.  A„  must  be  a  periodic  function  of  the 
wavenumber  k0  =  2ir/c0-  For  l«l  »  1,  |  A„  /  Ac  |  be¬ 
comes  quite  small,  as  shown  in  the  figure.  A  small  value 
of  A„  does  not  necessarily  mean  that  the  far  shifted  peri¬ 
odic  modes  are  not  important.  This  is  because  the  y"  terms 
in  (8)  get  large  negative  values  for  |  n  |  »  1 .  w  >  u>r 
and  u/k.  <  c.  Under  this  condition  of  imaginary  argu¬ 
ments.  the  Bessel  functions  70  are  replaced  by  the  modi¬ 
fied  Bessel  functions,  /0.  This  fact  results  in  large  value 
°f  A) ( I  yn  I )  in  (15)  when  we  superpose  electric  fields  of 
orthogonal  modes  at  the  boundary  of  the  slow  wave  struc¬ 
ture.  In  other  words.  Dm„A„  decrease  very  slowly  as  |  n  | 
increases,  because  |  v„ or  |  is  not  smaller  than  unity. 

For  example,  given  k.  =  4.01  cm"1  and  Np  =  2  x  10" 
cm"  ',  one  obtains; 

(D_4._4)M_4)  =  (1.5  x  10(’)  x  (8.0  x  1 0 “ 5 )  =  120 
(Dn.0)Mn)  =  (27)  x  (1.00)  =  27 
(Daa)(A4)  =  (5.7  x  10")  x  (1.3  x  10"10)  =  75 

(D.aa)(A4)  =  (-3.2  x  I09)  x  (1.3  x  10'"’) 

=  -0.42 

(D4._4)M_4)  =  (-2.7  x  10J)  x  (8.0  x  10'5) 

=  -0.21. 

This  means  that  the  relative  magnitude  of  the  off-diagonal 
terms  of  DmnA„  are  not  negligibly  small  in  comparison  to 
the  major  diagonal  terms.  The  right-hand  side  of  D  =  0 
is  of  the  order  of  10':  in  practical  computations.  In  con¬ 
clusion.  we  must  superpose  many  displaced  modes  to  sat¬ 
isfy  the  boundary  condition  (10)  when  a  -  I . 

B.  Dispersion  Relation  for  a  Solid  Beam  with  Finite 
Density 

So  far,  the  numerical  calculations  were  made  under  the 
assumption  that  Nh  =  0  and  the  growth  rate  Im(w)  =  0. 
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Nt«  2  «  to’^cm"5),  h»0.443(em) 


Fig.  5.  Eigenfrequencies  for  three  modes  versus  plasma  density  for  beam 
v/c  -  0.9. 


WAVENUMBER  k,  (em'1) 

Fig.  6.  Relative  magnitudes  of  orthogonal  modes.  versus  wave¬ 

number*.  for Nn  =  2  x  10"  cm-'. 

To  obtain  the  growth  rates  in  the  presence  of  the  beam, 
nonzero  values  of  ufc  must  be  included— we  are  now  look¬ 
ing  for  complex  roots  of  u  for  real  values  of  the  wave- 
number,  k.,  assuming  practical  values  of  the  beam  density 
Nb.  In  the  following  calculations,  a  solid  beam  with  the 
radius  Rh  =  0.9  cm  and  a  density  of  Nh  =  5.15  x  1010 
cn."3  was  assumed. 

The  calculation  was  performed  as  follov  we  plot  a 
mapping  of  absolute  values  j  D  |  of  the  9  x  9  determinant 
D,  on  the  i  plex  plane  for  given  parameters.  Because 
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An  example  of  the  contour  mapping  is  shown  in  Fig. 

7.  where  Np  =  2  x  10"  cm"1,  Nh  =  5.15  x  1010  cm"  , 
v/c  =  0.9,  and  k.  =  2.0  cm'1  have  been  assumed.  In 
this  figure,  the  horizontal  and  vertical  axes  are  Re(u)/2* 
and  lm(u)/2ir,  respectively.  Three  real  roots  and  a  pair 
of  complex  conjugate  roots  are  observed  at  the  centers  of 
circles.  The  heart-shaped  white  part  near  in  the  middle  of 
the  figure  is  where  the  values  of  the  determinant  are  ex- 
tremelv  large  (for  example,  on  the  order  of  1040).  The 
disper  >n  relation  for  complex  u  versus  r  ,al  k-  for  Nn  = 

2  x  10“  cm"3,  Nh  =  5.15  x  10u'  anc  v/c  =  0.9  is 
shown  in  Fig.  8.  The  thick  solid  li  -  are  the  real  part  of 
the  complex  freo  cy  u/2x,  and  the  thin  solid  lines  are 
for  the  real  root  The  result  is  not  very  different  from 
those  in  [13]  and  [14]  for  vacuum  BWO’s,  except  that  the 
growth  rates  are  a  function  of  plasma  density  in  this  case. 

The  two  thin  lines  close  to  v/c  =  0.9  are  the  slow 
(SBSCW)  and  fast  (FBSCW)  beam  space-charge  waves, 
as  was  shown  in  [14]. 

The  dispersion  relation  for  relatively  large  plasma  den¬ 
sity,  Np  =  2  x  1012  cm"3,  and  t\'h  =  5.15  x  1010  cm"3 
is  shown  in  Fig.  9.  Interactions  between  BSCW  and  the 
A  or  B  modes  are  depicted  in  detail.  There  are  again  the 
SBSCW  and  the  FBSCW  near  the  v/c  =  0.9  line.  The 
peak  growth  rate  Imfta)^  can  be  found  near  the  point 
where  the  complex  conjugate  roots  and  two  real  roots 
merge. 

The  imaginary  part  Im(u)/2x  versus  k.  is  plotted  for 
various  conditions  in  Fig.  10.  The  solid,  chained,  and 
dashed  lines  correspond  to  the  intersections  of  the  BSCW 
with  the  TM0(,  A,  and  B  mode  respectively.  The  cases 
I,  II,  and  III  are  for  Np  =  2  x  i  2  x  10l5,  and  3.2  x 
1012  cm"3,  respectively.  For  relatively  small  values  of  Np 
of  the  case  I,  the  dominant  mode  of  oscillation  is  the  TM0, 
mode,  whereas  for  large  values  of  Np,  i.e.,  cases  II  and 
III,  the  growth  rates  for  B-mode  approach  to  those  for 
TM0j  mode,  resuitii  .n  a  possibility  of  higher  frequency 
mode  to  oscillate. 

In  Fig.  11  the  peak  growth  rate,  Im(u))ma,/2T  is  plot¬ 
ted  versus  the  plasm  density  Np  in  the  slow  wave  struc¬ 
ture.  The  growth  rates  do  not  change  appreciably  up  to 
densities  of  Np  =  6  x  1010  cm"3.  In  such  cases,  the  max¬ 
imum  growth  rate  for  the  TM0I  mode  is  considerably 
greater  than  those  of  the  A  and  B  modes.  Around  Np  — 

1012  cm"3,  the  pet  srowth  rate  for  the  B  mode  ap¬ 
proaches  those  for  ti  .  M0i  mode.  In  other  words,  other 
high-frequency  oscillations  may  arise  when  V„  is  high,  in  **■ 
addition  to  TM0i  mode  with  an  increased  oscillation  fre¬ 
quency.  This  is  main  result  obtained  in  the  present 
linear  analysis. 


£>(w,  k.)  s  D( u0,  k.)  +  (w  -  uo)  —  =  K  (16) 

flu 

and  D(  u0,  k.  )  =0,  the  solution  can  be  found  at  the  center 
of  the  contour  circle  lines  of  j  D  |.  This  is  because  (16)  is 
an  equation  of  circles  with  radius  |  AT|  in  the  complex  u 
plane. 


IV.  Discussion  and  Conclusion 
We  have  numerically  analyzed  the  linear  dispersion  re¬ 
lation  of  a  plasma-loaded  corrugated-wall  slow  wave 
structure  immersed  in  an  infinitely  large  axial  magnetic 
field.  It  is  found  that  the  terms  including  L„  in  (13)  affect 
only  Im(u);  Re(u)  is  unchanged  even  if  the  £*  term 
neglected. 
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Fig.  7.  Contour  map  of  absolute  value  of  9  x  9  determinant  |D|  on  complex  w/2*  plane  for  A/,,  =  2  x  10"  cm  .  N,,  -  5. 15 
x  10"’  cm  h  -  0.445  cm.  c,,  =  1 .67  cm.  and  k.  =  2.0  cm'1.  Three  real  roots  and  pair  of  complex  conjugate  roots  are 
shown. 


Np=  2  x  1o"(cm"3)  .  Nb=  5.15  x  1o'°(em’3) 


Fig.  8.  Dispersion  curves  and  linear  growth  rates  for  Np  =  2  x  I0"cm  \ 
N„  «  5.15  x  I0‘° cm'1. 


Fig.  10.  Growth  rale.  Im(w)/2T.  versus  wavenumber*..  Case  I:  Nr  =  2 
x  10"  cm  ':  Case  II:  N,  -  2  x  10’’  cm  ':  Case  111:  N„  =  3.2  x  10'-’ 
cm'1.  In  all  cases.  Nh  =  5.15  x  10"’  cm  ’. 


Fig.  9.  Dispersion  curves  for  N,,  =  2  x  I0';  cm  '.  Nh  =  5.15  x  10*" 
cm  '.  Only  real  parts  are  shown  in  detail. 


Nfc=  5.15x1  ty°{  c  m' *  ) 


PLASMA  DENSITY  Np{*10  cml 

Fig.  I).  Peak  growth  rates  Im («)„,, /2»  versus  plasma  density. 

It  was  found  that  for  deep  corrugation  (a  =  0.308), 
corresponding  to  our  practical  BWO  devices,  nine  dis¬ 
placed  periodic  modes  must  be  superposed,  and  that  terms 
up  to  the  order  of  (yBcr)10  in  (14)  must  be  calculated  to 
have  numerical  results  within  an  accuracy  of  1  % . 

According  to  the  linear  theory,  the  primary  effect  of  the 
presence  of  plasma  in  the  slow  wave  structure  is  an  in- 
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crease  in  the  oscillation  frequency  of  the  TM0|  mode.  The 
second  effect  is  an  appearance  of  higher  modes.  For 
plasma  densities  greater  than  1012  cm'3,  an  oscillation  of 
B  mode  can  be  predicted,  as  s.’iown  in  Fig.  11.  A  mode¬ 
switching  observed  in  our  experiment  [8]  may  be  related 
to  the  result  of  the  present  linear  theory. 

The  enhancements  of  the  oscillation  observed  in  our  ex¬ 
periments  may  be  h  -  ond  the  scope  of  the  linear  theory 
derived  in  this  papei  The  experimental  observations  in¬ 
dicate  a  strong  enhancement  at  the  TM0)  interaction.  It  is 
possible  that  a  nonlinear  process  is  responsible  fc :  •  ob¬ 
served  enhancement.  Further  a'  ’lyses  to  undersu-.u  the 
observed  enhancement  are  undt  .vay. 
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A  novel  method  for  calculating  the  dispersion  relation  of  electromagnetic  modes  in  an 
arbitrary  periodic  slow  wave  structure  is  reported.  In  this  method  it  is  sufficient  to  know  the 
frequencies  corresponding  to  three  special  wave  number  values,  with  other  points 
calculated  using  an  approximate  analytical  expression.  This  technique  was  successfully 
applied  to  determine  the  dispersion  relation  of  the  TM^  mode  in  a  sinusoidally  corrugated 
slow  wave  structure.  This  structure  is  commonly  used  in  relativistic  high-power  backward 
wave  oscillators  and  traveling-wave  tubes,  and  is  expected  to  have  many  additional 
applications. 


When  an  electron  beam  is  injected  into  a  periodic  slow 
wave  structure,  the  beam  structure  resonance  often  leads  to 
either  an  absolute  or  a  convective  instability.  The  beam 
space-charge  waves  couple  to  the  slow  wave  structure  nor¬ 
mal  modes  to  produce  microwave  and  millimeter  wave  ra¬ 
diation.  The  group  and  phase  velocity  of  the  slow  wave 
structure  modes,  as  well  as  the  electron  beam  characteris¬ 
tics,  determine  the  nature  and  frequency  of  the  beam-wave 
interaction. 

There  are  many  families  of  microwave  and  millimeter 
wave  generating  devices  whose  operation  depends  on  this 
type  of  interaction.  As  an  example.  Fig.  1  shows  schemat¬ 
ically  the  regions  of  operation  of  some  of  these  devices  in 
the  frequency-wave  number  domain,  as  well  as  some  typi¬ 
cal  slow  wave  structures.  Relativistic  traveling-wave  tubes1 
(RTWTs)  operate  below  the  point  where  the  wave  number 
(normalized  to  the  structure  period)  is  equal  to  v  radians. 
In  this  case,  both  the  phase  and  the  group  velocities  are 
positive.  At  PL  —  tr,  the  upper  cutoff  frequency,  the  elec¬ 
tromagnetic  wave  undergoes  a  phase  shift  of  v  radians  per 
period  of  the  slow  wave  structure.  At  this  point,  the  phase 
velocity  is  positive  and  the  group  velocity  is  zero.  Back¬ 
ward  wave  oscillators  (BWOs)  as  well  as  cardnotrons  op¬ 
erate  in  the  region  tr <pL <2v,  where  the  phase  velocity 
(ai/P)  is  positive  and  the  group  velocity  (do) /dp)  is  neg¬ 
ative  (<u  is  the  angular  frequency  and  P  is  the  wave  num¬ 
ber).  Extended  interaction  oscillators  (EIO)  operate  very 
close  to  the  PL  =  2ir  point. 

Relativistic  backward  wave  osillators2  and  related  de¬ 
vices  which  typically  operate  close  to  the  PL  =  tr  point, 
prove  to  be  efficient  and  powerful  microwave  and  millime¬ 
ter  wave  sources  with  reported  record  power  levels  reach¬ 
ing  IS  GW  at  wavelength  of  3  cm3  and  S  GW  at  3  mm,4 
with  efficiencies  of  up  to  50%.  In  these  devices,  a  smoothly 
corrugated  waveguide  has  been  found  to  be  advantageous 
over  alternative  periodic  structures  due  to  its  high-power 
handling  capabilities.  Many  other  devices  (to  be  discussed 
later)  also  employ  corrugated  waveguide  slow  wave  struc¬ 
tures  operating  in  both  TM  and  TE  modes.  It  is  important, 
therefore,  to  know  the  dispersion  relation  of  such  slow 
wave  structures  in  order  to  synchronize  the  phase  velocity 
of  the  wave  with  the  electron  beam  to  produce  efficient 


interaction.  In  this  letter  we  present  a  novel,  accurate,  and 
simple  method  for  calculating  the  dispersion  curve  of  an 
infinitely  long  slow  wave  structure  of  arbitrary  geometry. 
This  method  was  successfully  applied  to  determine  the  dis¬ 
persion  curve  of  the  TM^,  mode  in  a  sinusoidally  corru¬ 
gated  waveguide.  The  results  are  then  compared  with  those 
achieved  by  using  other  techniques. 

Any  periodic  slow  wave  structure  with  n  periods,  when 
shorted  at  both  ends,  will  exhibit  (n  f  1)  resonant  fre¬ 
quencies  with  a  phase  shift  per  period  equally  spaced  be¬ 
tween  0  and  v.  Outside  of  the  region  0</?L<ir,  the  disper¬ 
sion  relation  is  periodic  in  P  space  (Floquet’s  theorem). 
This  dispersion  relation  can  be  either  be  calculated  using 
computational  methods  or  measured  experimentally  (cold 
test). 


(c) 


FIG.  1.  (a)  Repots  of  operation  in  the  frequency-wavelength  domain  of 
various  microwave  and  millimeter  wave  generators;  (b)  schematic  dia¬ 
gram  of  a  sinuaoidally  corrugated  waveguide;  (c)  iris  loaded  waveguide. 
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TABLE  1.  Resonance  frequencies  (GHz)  J  a  sinusoidally  corrugated 
waveguide  cavity  ( R0  =  1.5  cm,  L  =  1.67  cm,  h  =  0.273  cm). 


Resonance  frequency 
(GHz) 

2-period  cavity 

4-period  cavity 

/,=/10*-  =  0) 

7.404 

7.412 

A  0L  =  */*) 

N/A 

7.589 

/„;=/TBL  =  r/2) 

8.046 

8.046 

A  &L-  =  3<r/4) 

N/A 

8.557 

fr=  A  0L  =  IT 

8.773 

8.754 

Instead  of  using  a  large  number  of  periods  to  accu¬ 
rately  calculate  the  dispersion  relation,  we  will  show  that  it 
is  sufficient  to  calculate  the  short  circuit  cavity  resonant 
frequencies  corresponding  to  three  special  normalized 
wave  number  values,  PL  =  0,  ir/2,  and  ir  radians  by  rigor¬ 
ously  solving  the  boundary  value  problem.  With  these 
three  frequencies,  one  can  use  an  analytical  expression  to 
produce  the  complete  dispersion  curve.  Furthermore,  we 
can  also  analytically  calculate  the  phase  and  group  veloc¬ 
ities  at  any  point  along  the  dispersion  relation.  We  have 
used  a  two-dimensional  electromagnetic  code5  ("Super¬ 
fish”  )  for  the  calculation  of  these  three  special  frequencies 
for  a  sinusoidally  corrugated  waveguide,  oriented  in  the  z 
direction  whose  radius  is  described  by 

/?=/?oll  +  h  cos(2irz/£)],  (1) 


FIG.  2.  (a)  Mode  patterns  (electric  field  lines)  corresponding  to  zero 
phase  shift,  (b)  ir/2  phase  shift,  and  (c)  v  phase  shift  in  the  TMoi  lower 
pass  band  of  a  sinusoidally  corrugated  waveguide. 


where  R0  is  the  average  radius,  h  is  the  normalized  corru¬ 
gation  amplitude,  and  L  is  the  structure  period. 

For  high-power  relativistic  BWOs,  the  main  mode  of 
interest  is  the  cylindrically  symmetric  TMo,  mode.  In  our 
test  case,  the  waveguide  average  radius  was  R0  =  1.50  cm, 
the  ripple  period  was  L=  1.67  cm,  and  the  normalized 
corrugation  amplitude  was  h  =  0.273  cm.  Superfish  mod¬ 
eling  yielded  (with  proper  boundary  conditions)  the  three 
special  resonance  frequencies  in  the  lower  passband.  The 
results  are  tabulated  in  the  second  column  of  Table  I  (two 
period  cavity).  The  corresponding  mode  patterns  (electric 
field  lines)  are  shown  in  Figs.  2(a),  2(b),  and  2(c),  which 
pertain  to  phase  shifts  per  period  of  0,  ir/2,  and  v,  respec¬ 
tively. 

Next,  the  expression  used  to  calculate  the  complete 
dispersion  relation  is  discussed.  By  using  the  impedance  or 
ABCD  matrix  of  a  four-terminal  network,  together  with 
Floquet’s  theorem6  to  describe  a  periodic  slow  wave  struc¬ 
ture,  the  dispersion  relation  car  be  expressed  in  the  non¬ 
explicit  form 

f =f(cos  pL,G)  -  (2) 


The  three  terms  in  Eq.  (3)  are  sufficient  to  ensure  excellent 
accuracy  for  most  practical  cases,  as  will  be  shown  later. 
Notice  that  the  geometrical  factor  G  (as  yet  unspecified) 
does  not  appear  in  the  equation.  Rather,  it  enters  indirectly 
through  A,  B,  C,  which  are  geometry-dependent  parame¬ 
ters  with  simple  physical  meaning.  The  first  term.  A,  rep¬ 
resents  the  value  of  the  frequency  near  the  midband.  The 
second  is  an  increment  term  whose  maximum  value  equals 
half  the  difference  between  the  upper  and  lower  cutoff  fre¬ 
quencies  of  the  passband.  The  third  term  is  a  correction 
term.  The  combination  of  all  three  terms  exactly  satisfied 
the  dispersion  relation  at  the  three  special  wave  numbers 
(0,  7r/2,  and  n)  and  is  an  excellent  approximation  at  all 
other  wave  numbers,  as  will  be  shown  later. 

Equation  (3)  enables  one  to  directly  calculate  the 
phase  ard  group  velocities  at  any  point  along  the  disper¬ 
sion  re  in.  The  expressions  derived  for  the  velocities  are 
given  in  £qs.  (4)  and  (5): 

[ILrrf  +  Wfnn)) 

fptuse  cos _ '(  [(B2-4CA/)l/i-R]/2Q’  (4) 


where  /is  the  frequency  of  the  ctromagnetic  .  .iation 
and  G  is  a  geometric  factor  rela;  to  a  specific  slow  wave 
structure  dimension.  Let  A  —f„n,  B  =  (f„—f0)/2, 
C = A/2  —  ( A  +  /o)/2.  It  can  be  shown  that  the  dispersion 
relation  of  any  general  periodic  structure  can  be  expressed 
in  the  following  explicit,  approximate  form: 

f =A—B  cos  (3L  —  C  cos2  PL.  (3) 


vgroUp=do}/dP—2irL(B  sin  PL  +  2C  cos  PL  sin  PL),  (5) 
where  A/ =f-fn/2. 

The  results  of  our  test  case  are  plotted  in  Fig.  3.  It 
shows  the  TM^  dispersion  curve  in  the  lower  passband  for 
the  sinusoidally  corrugated  waveguide  Fig.  3(a),  as  well  as 
that  of  an  iris  loaded  waveguide  having  the  same  minimum 
and  maximum  radial  dimensions.  The  two  are  similar  with 
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FIG.  3.  (a)  TM,ii  dispersion  curve  for  a  sinusoidal!,  corrugated,  iris 
loaded,  and  smooth  waveguide,  (b)  normalized  group  velocities,  (c) 
phase  velocities. 


the  iris  loaded  guide  dispersion  curve  always  lower  in  fre¬ 
quency  compared  to  the  corrugated  guide.  This  is  expected 
since  the  volume  of  the  first  is  larger  than  that  of  the 
second,  driving  all  the  resonant  frequencies  down.  The 
open  circles  at  (JL  =  0,  tr/2,  and  ir  were  calculated  by  using 
Superfish  and  were  used  in  Eq.  (3)  to  calculate  the  com¬ 
plete  dispersion  relation.  For  comparison,  the  dispersion 
relatioA  of  a  smooth  waveguide  is  also  given.  Figures  3(b) 
and  3(c)  display  the  phase  and  group  velocities  (normal¬ 
ized  to  the  speed  of  light)  for  the  same  cases  [Eqs.  (4)  and 
(5)]. 

In  order  to  estimate  the  accuracy  of  this  technique  as 
applied  to  the  sinusoidally  corrugated  guide,  we  first  cal¬ 
culate  additional  resonance  points  on  the  dispersion  curve, 
using  “Superfish'*.  Four  periods  of  the  corrugated  guide 
cavity  were  modeled,  yielding  five  resonances  in  the  wave 
number  range  of  0  to  ir  and  the  results  are  given  in  the 
third  column  of  Table  I  (four-period  cavity)  and  are  plot¬ 
ted  as  additional  open  circles  in  Fig.  3(a)  with  very  good 
agreement  (eight  periods  were  also  successfully  modeled). 
As  an  independent  check,  the  dispersion  relation  derived 
here  was  compared  to  the  dispersion  relation  deived 
elsewhere78  for  the  same  corrugated  guide  geometry  by 
solving  a  system  of  coupled  linear  differential  equations. 
Figure  4  shows  a  companion  between  the  dispersion  curve 


9.or 


FIG.  4.  Comparison  of  the  TMoi  dispersion  curve  of  Ref.  7  (solid  line) 
and  this  work  (dotted  line)  for  a  sinusoidally  corrugated  waveguide. 

as  calculated  by  these  two  techniques  [Ref.  7  and  Eq.  (3)] 
for  the  same  corrugated  guide.  The  results  are  in  excellent 
agreement  (within  0.15%)  over  the  entire  wave  number 
range  0  <(3L<it.  Outside  this  range  the  dispersion  relation 
is,  of  course,  periodic. 

It  is  expected  that  this  technique  will  be  beneficial  for 
determining  the  dispersion  curve  of  periodic  structures  of 
arbitrary  shapes.  The  three  special  resonance  frequencies 
which  are  needed  in  order  to  generate  the  complete,  gen¬ 
eral  dispersion  curve  can  be  calculated  by  rigorously  solv¬ 
ing  the  boundary  problem  or  measured  experimentally. 

Corrugated  slow  wave  structures  may  prove  beneficial 
in  prospective  applications  such  as  TE  mode  slow  wave 
cyclotron  amplifiers,9  TM  mode  relativistic  extended  inter¬ 
action  oscillators  and  amplifiers,10  TE  mode  CARMs,11 
electromagnetically  pumped  free-electron  lasers,1213  and 
plasma-loaded  backward  wave  oscillators.14 

This  work  was  sponsored  in  part  by  the  Weapons  Lab¬ 
oratory  (Kirtland  AFB)  through  a  contact  administered 
by  NRL. 
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An  experimental  demonstration  of  a  strong  enhancement  of  the  interaction  efficiency  in  a  high-power 
relativistic  backward-wave  oscillator  when  plasma  is  injected  is  presented.  Controlled  plasma  injection 
enhances  the  interaction  efficiency  over  the  vacuum  case  by  a  factor  of  up  to  8  to  a  value  of  about  40%. 
The  enhanced  interaction  is  attributed  to  induced  scattering  of  the  electromagnetic  radiation  of  electrons 
in  an  electrostatic  field  produced  by  the  background  plasma  and  a  beam  space-charge  wave  in  the  corru¬ 
gated  interaction  region. 

PACS  numbers:  52.60.+h,  42.55.Tb,  52.40.Mj.  52.70.Gw 


Many  varieties  of  high-power  microwave  (HPM)  de¬ 
vices  have  been  studied  in  recent  years,  including  the 
well-known  magnetron  and  klystron  as  well  as  newer  de¬ 
vices  such  as  the  gyrotron,  free-elcctron  laser,  and  virtual 
cathode  oscillator. 1  All  of  these  sources  have  one  thing 
in  common— they  are  driven  by  an  intense,  monoener- 
getic,  unneutralized  electron  beam  which  unstably  in¬ 
teracts  in  high  vacuum  with  an  electromagnetic  wave, 
leading  to  the  conversion  of  the  beam’s  kinetic  energy 
into  electromagnetic  radiation.  The  power  levels  avail¬ 
able  from  such  devices  have  grown  by  an  order  of  magni¬ 
tude  every  decade  since  1940,  reaching  about  10'°  W  at 
a  wavelength  of  3  cm  arm  10*  W  at  3  mm.2 

This  work  demonstrates  for  the  first  time  a  strong 
efficiency  enhancement  in  a  relativistic  backward-wave 
oscillator  (BWO)  by  external  plasma  injection.  Plasma 
effects  in  conventional  microwave  devices  can  usually  be 
neglected  since  the  plasma  frequency  of  the  background 
gas  is  much  smaller  than  the  plasma  frequency  of  the 
electron  beam  <  10  ~2).  Recent  theoretical 

studies,1  however,  have  predicted  that  the  presence  of  a 
plasma  in  HPM  devices  may  lead  to  greatly  enhanced 
performance,  attracting  renewed  scientific  interest4  The 
presence  of  a  background  plasma  can  serve  to  increase 
the  space-charge  limiting  current  by  a  factor  of 
(1  where  f^ni/n,  represents  any  neutralization 

provided  by  the  ions  in  the  background  plasma.  The 
space-charge  limiting  current  in  the  presence  of  the  plas¬ 
ma  for  a  thin  hollow  beam  with  a  mean  radius  a  and  a 
thickness  A«a  is  given  by 


nlr^-D* 
"  l21n(6/a)JIl  -/I 


(1) 


Here  yp  is  the  relativistic  mass  factor  for  the  beam  elec¬ 
trons,  and  b  is  the  radius  of  the  drift  tube.  By  increasing 
the  beam  current  through  introduction  of  a  background 
plasma,  a  new  limitation  is  encountered  which  is  imposed 
by  the  onset  of  beam-plasma  instability.  This  instability 
will  occur  at  current  levels  which  are  larger  than  the  vac¬ 


uum  limiting  current  by  the  factor1 

rl(i-r'2)/Ci-r'2/,)l3/2.  (2) 

In  the  case  of  mildly  relativistic  electron  beams  having 
y“3,  the  actual  current  carried  by  the  beam  in  the  pres¬ 
ence  of  plasma  may  be  7  times  larger  than  the  vacuum 
case.  The  increased  injected  current  level  simply  allows 
operation  at  higher  beam  power  and  associated  higher 
microwave  power  output  without  affecting  the  interac¬ 
tion  efficiency  or  the  physical  interaction  mechanism. 

The  plasma  presence  in  the  device  can,  however,  also 
completely  alter  the  nature  of  the  interaction  mecha¬ 
nism,  leading  to  greatly  enhanced  device  efficiency.  The 
studies  reported  here  are  the  first  which  clearly  belong  to 
this  latter  category. 

Since  BWO’s1  are  simple  devices,  provide  fairly 
effective  conversion  of  electron  beam  energy  into  radia¬ 
tion,  and  can  easily  be  filled  with  plasma,  they  are  ideal 
candidates  for  the  evaluation  of  plasma  effects.  Initial 
studies6,7  utilized  background  plasmas  produced  by 
electron-beam-impact  ionization  of  a  low-pressure  neu¬ 
tral-gas  background,  and  demonstrated  enhanced  mi¬ 
crowave  output  powers.  Interpretation  of  the  results, 
however,  and  identification  of  the  physical  mechanisms 
involved  have  proven  difficult  In  some  of  these  earlier 
studies,6  the  increased  output  power  was  attributed  to  a 
substantial  increase  in  injected  beam  current  allowed  by 
the  increase  in  the  space-charge  limit  allowed  by  the 
plasma.  In  another,7  beam  current  was  again  increased 
over  the  vacuum  BWO  studies,  and  identification  of  the 
mechanism  was  complicated  by  the  strong  resonant 
dependence  of  microwave  output  oo  the  externally  ap¬ 
plied  magnetic  field.  In  the  work  presented  here,  inject¬ 
ed  beam  parameters  were  carefully  held  constant  for 
both  vacuum  and  plasma-filled  BWO  experiments.  The 
resulting  enhancement  in  BWO  microwave  output  power 
when  plasma  is  present  is  therefore  dearly  a  result  of  an 
improvement  in  the  electronic  efficiency  of  the  device. 
Thus  this  work  unambiguously  demonstrates  efficiency 
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FIG.  I.  Schematic  diagram  of  a  relativistic  backward-wave  oscillator  with  external  plasma  injection. 


enhancement  in  a  high-power  microwave  tube. 

In  the  experiments,  an  independently  controllable  ar¬ 
gon  plasma  source  was  used  to  externally  inject  plasma 
directly  into  the  BWO  structure  as  shown  in  Fig.  1.  A 
hollow,  relativistic  electron  beam  (630  kV,  13  kA,  100 
ns)  of  average  radius  0.8  cm  and  electron  density  of 
about  5xio"  cm-3  was  injected  into  a  BWO 
corrugated-wall  slow-wave  structure  immersed  in  a  uni¬ 
form  axial  magnetic  field  of  about.  12  kG.  A  coaxial 
p'asma  gun,'  located  about  100  cm  downstream  in  a 
field-free  region,  generated  an  argon  plasmoid  which 
crossed  magnetic  field  lines  at  an  average  velocity  of 
about  1.2  cm//ii  on  its  way  towards  the  interaction  re¬ 
gion.  The  entire  system  was  evacuated  to  pressure 
<  4x  10  ~5  Ton.  The  plasma  column  parameters  were 
measured  with  35-  and  70-GHz  microwave  interferome¬ 
ters  (for  density)  and  Langmuir  probes  (for  velocity, 
density,  and  temperature).  Initial  results  indicate  that 
the  average  plasma  density  could  be  varied  continuously 
by  changing  the  gun  voltage  and  gas  pressure  from  zero 
up  to  a  maximum  value  estimated  at  ~1013  cm-3, 
which  is  well  above  the  corresponding  beam  density 
(5xio"  cm-3). 

Efficiency  enhancement  was  observed  over  a  wide 
range  of  injected  plasma  densities,  0  <  n  <  na,  where  na 


FIG.  1  Peak  power  output  at  8.4  GHz  vs  the  time  interval 
between  the  plasma  gun  operation  and  the  electron  beam 
firing. 


is  a  critical  plasma  density.  This  critical  plasma  density, 
to  be  discussed  later,  occurred  at  Ta<90  /is  on  the  plas¬ 
ma  pulse  time  scale.  The  efficiency  enhancement  factor 
over  vacuum  BWO  operation  was  found  to  be  density 
dependent,  and  reached  a  maximum  value  of  8  when  the 
electron  beam  was  injected  into  an  optimized  plasma 
density.  This  occurs  twice  on  the  plasma  pulse  time 
scale,  once  during  density  rise  (7*a« 60  /is)  and  then 
again  during  density  decay  (Tat  100  /is);  see  Fig.  2.  At 
these  points  the  interaction  efficiency  increased  to  almost 
40%  com oa red  with  about  5%  for  the  vacuum  BWO  un¬ 
der  the  s  operating  conditions. 

From  c  considerations  it  is  possible  to  estimate  the 
critical  piamia  density  which  will  overdrive  the  device 
and  ouench  the  interaction.  Figure  3  depicts  the  approx¬ 
imate  di  rsion  relation  of  the  lowest -order  symmetric 
electromagnetic  branches9  (labeled  TMoi  and  TM02)  as 
well  as  the  plasma  branenes  (labeled  plasma  waves) 
which  we  shall  refer  to  later.  Note  that  the  intersection 
point  of  the  TMoi  mode  and  the  electron-beam  space- 
charge  wave  lies  on  a  portion  of  the  electromagnetic- 
wave  dispersion  relation  having  negative  group  velocity. 


FIG.  3.  The  dispersion  relation  of  plasma-filled  corrugated 
waveguide  including  the  electromagnetic  and  plasma  branches, 
with  superimposed  ideal  beam  wave  added.  K»  is  the  beam  ve¬ 
locity,  e  is  the  speed  of  tight,  and  6*1.445  cm  is  the  average 
radius  of  the  corrugated -waveguide,  with  radius  given  by 
r  *6+0.445sin(2*r/l.67)  cm. 
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which  implies  a  backward-propagating  wave.  The 
backward-propagating  wave  is  reflected  by  a  waveguide 
beyond  cutoff  at  the  injection  point  and  then  detected 
downstream  from  the  structure  (Fig.  1).  This  intersec¬ 
tion  occurs  close  to  the  upper  cutoff  frequency  of  the 
TMoi  passband  at  8.4  GHz  (point  A  in  Fig.  3).  To  first 
order,  the  presence  of  a  magnetized  plasma  in  the  slow- 
wave  structure  will  raise  its  lower  cutoff  frequency  above 
the  vacuum  value  of  7.3  GHz  (point  B).  If  the  device  is 
overdriven  to  the  point  where  the  lower  cutoff  frequency 
in  the  presence  of  plasma  is  equal  to  or  greater  than  the 
upper  cutoff  frequency  in  the  absence  of  plasma,  no 
backward-wave  interaction  is  possible. 

An  approximate  expression  for  the  cutoff  frequency 
(TM  modes)  of  a  waveguide  filled  with  magnetized  plas- 
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In  Eq.  (3)  oip  is  the  background  plasma  frequency,  c  is 
the  speed  of  light,  are  the  roots  of  the  Bessel  func¬ 
tions  corresponding  to  the  mode  of  interest,  and  b  is  the 
waveguide  radius.  When  solved  for  the  critical  back¬ 
ground  plasma  density  which  would  quench  the  interac¬ 
tion,  using  our  experimental  parameters,  Eq.  (3)  yields 


/ICT-2X10"  cm-3. 


(4) 


When  the  backward-wave  oscillator  was  overdriven 
with  plasma  above  the  critical  density,  the  fundamental 
TMoi  mode  at  8.4  GHz  was  quenched  and  strong, 
higher-frequency  emission  (12-18  GHz)  was  observed, 
possibly  indicating  mode  switching  from  TMoi  to  TM02 
(point  C  in  Fig.  3). 

The  dramatic  mode  switching  is  shown  in  Fig.  4.  It 
shows  the  amplitude  (vertical  scale)  of  different  frequen¬ 
cy  (horizontal  scale)  components  of  the  backward-wave 
oscillator  using  the  dispersive  line  technique  (different 
frequencies  are  resolved  by  their  different  propagation 
velocities  in  a  dispersive  waveguide).  Figure  4(a)  repre¬ 
sents  the  spectrum  emitted  by  the  backward-wave  oscil¬ 
lator  under  the  condition  of  strong  plasma  enhancement, 
where  the  dominant  frequency  component  is  near  8.4 
GHz.  An  undispersed,  prompt  microwave  output  signal 
is  superimposed  on  the  dispersed  signal  as  a  timing  refer¬ 
ence.  Figure  4(b)  represents  the  spectrum  for  an  over¬ 
driven  device,  in  which  case  mode  switching  occurred 
and  frequency  components  in  the  range  10-18  GHz  are 
present  (TM02  or  other  high-order  modes  as  well  as  plas¬ 
ma  oscillation).  At  even  larger  plasma  densities,  mi¬ 
crowave  breakdown  within  the  device  will  occur. 

Strong  microwave  emission  was  also  detected  in  the 
18-26-GHz  band,  proportional  in  amplitude  to  the  fun¬ 
damental  TMoi  backward-wave  oscillations.  These  high 
frequencies  are  believed  to  be  produced  by  a  free- 
electron  laser  interaction  driven  by  an  electromagnetic 
pump,  as  reported  previously.12'11 
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FIG.  4.  Spectral  results  from  plums -enhanced  backward- 
wave  oscillator,  (a)  Maximum  enhancement  and  (b)  over¬ 
driven  device. 

As  these  experiments  demonstrate,  the  controlled  pres¬ 
ence  of  a  plasma  in  the  corrugated  structure  strongly 
enhances  the  interaction  efficiency  over  a  very  wide 
range  of  plasma  densities.  To  understand  the  physical 
mechanism  responsible  for  the  enhanced  interaction.  Fig. 
3  also  depicts  the  approximate  dispersion  relation  of  a 
plasma-filled  corrugated  waveguide.  It  consists  not  only 
of  the  known  slow  electromagnetic  branches,  but  the 
plasma  branches  as  well  (labeled  plasma  waves  in  Fig. 
3).  These  low-frequency  modes  are  purely  plasma  waves 
and  vanish  in  the  limit  of  small  background  plasma  den¬ 
sity  (<dp— 0).  They  are  the  corrugated-waveguide  ver¬ 
sion  of  the  T rivelpicce-Gould  (TG)  modes."  Normally, 
those  plasma  oscillations  are  trapped  in  the  plasma  and 
cannot  be  coupled  out  efficiently  in  the  form  of  elec¬ 
tromagnetic  radiation.  However,  for  the  range  of  values 
of  the  plasma  density  and  beam  energy  relevant  to  our 
work,  the  beam  is  simultaneously  in  synchronism  (and 
therefore  can  exchange  energy)  with  the  backward 
branch  of  electromagnetic  wave  (point  A)  and  the  back¬ 
ward  branch  of  the  plasma  wave  (point  D  in  Fig.  3). 
Under  these  conditions,  all  three  waves  involved  have  the 
same  phase  velocity,  a>!k,  and  a  strong  component  of  axi¬ 
al  electric  field,  Et.  Therefore,  a  strong  enhancement  is 
expected.  It  is  believed  that  the  greatly  enhanced 
efficiency  observed  results  from  this  mechanism,  namely, 
induced  scattering  of  the  electromagnetic  radiation  of 
electrons  in  an  electrostatic  field  produced  by  the  back¬ 
ground  plasma  in  the  presence  of  the  corrugated  wall. 

The  introduction  of  plasma  into  the  device  did  not 
affect  the  diode  voltage  or  current,  and  the  beam  current 
actually  entering  the  slow-wave  structure  was  unaffected 
by  the  plasma.  No  diode  shorting  due  to  the  injected 
plasma  was  observed,  and  the  performance  of  the  plasma 
BWO  with  respect  to  variations  in  the  applied  axial  mag¬ 
netic  field  was  similar  to  the  vacuum  BWO. 14  No  sharp 
resooance  (similar  to  the  one  in  Ref.  7)  is  expected  for 
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this  mechanism,  and  none  was  observed. 

An  efficient  interaction  also  implies  a  very  broad  gain 
curve  with  respect  to  the  electron  beam  energy.  This  is 
expected  for  the  present  mechanism  since  the  electrons 
tend  to  stay  in  synchronism  with  both  backward 
branches  (points  A  and  D  in  Fig.  3)  even  as  the  beam  is 
losing  energy.  Although  the  effect  of  the  plasma  on 
beam  transport  inside  the  slow-wave  structure  could  not 
be  directly  measured,  no  evidence  of  significant  beam 
current  striking  the  structure  wails  was  observed. 

It  is  anticipated  that  plasma  injection  may  also  prove 
beneficial  in  other  high-power  micro*  <tve  devices  such  as 
free -electron  laser  and  gyrotrons,  prc  ding  frequency  tu- 
nability,15  overcoming  space-charge  limitations, 16  and 
enhancing  interaction  efficiencies. 17 
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